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ABSTRACT 


For  Terahertz  (THz)  imaging  to  be  useful  outside  of  a  laboratory  setting,  inexpensive,  yet 
sensitive,  detectors  such  as  uncooled  microbolometers  will  be  required.  Metamaterials 
can  improve  THz  absorption  without  significantly  increasing  the  thermal  mass  or  using 
exotic  materials  because  their  absorption  is  primarily  dependent  on  the  geometry  of  the 
materials  and  not  their  individual  optical  properties.  Finite  Element  (FE)  simulations 
revealed  that  an  array  of  squares  above  a  ground  plane  separated  by  a  dielectric  is 
efficient,  yet  thin.  Metamaterials  were  fabricated  and  their  absorption  characteristics  were 
measured  using  a  Fourier  Transfonn  Infrared  Spectrometer  indicating  that  the  FE 
simulations  are  accurate.  Metamaterial  structures  tuned  to  a  quantum  cascade  laser  (QCL) 
illuminator  were  incorporated  into  a  bi-material  sensor,  which  was  used  for  detection  of 
THz  radiation  from  the  QCL  source  with  good  sensitivity.  In  the  case  of 
microbolometers,  a  bolometric  layer  needs  to  be  embedded  in  the  metamaterial  to  form  a 
thin  microbridge.  Simulations  indicated  that  if  the  bolometric  layer  was  resistive  enough 
or  close  enough  to  the  ground  plane,  then  absorption  would  be  largely  unaltered. 
Metamaterials  with  a  conductive  Titanium  layer  embedded  into  the  dielectric  spacer  were 
fabricated  and  measured  with  an  FTIR,  confirming  this  behavior. 
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temperature  coefficient  of  resistance 
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I.  INTRODUCTION 

A.  TERAHERTZ  IMAGING 

In  recent  years,  there  has  been  a  surge  of  renewed  interest  in  the  terahertz  (THz), 
or  Far-Infrared  (Far-IR),  spectral  range.  The  THz  band  (Figure  1),  so  named  because  the 
frequencies  involved  are  typically  between  0.1-10  x  1012  Hz  (3  mtn-30  pm  wavelength), 
has  some  interesting  real  world  imaging  applications.  Terahertz  easily  penetrates  dry, 
non-polar,  low-conductivity  materials  (such  as  common  packaging  materials  and 
clothing),  but  is  readily  reflected  by  metallic  objects  and  strongly  absorbed  by  water.  This 
limits  the  effective  range,  especially  above  1  THz,  to  short-range  applications. 
Additionally,  THz  radiation  is  non-ionizing  and  does  not  penetrate  deeply  into  the  skin, 
making  it  potentially  safer  than  X-ray  imaging  technologies  and  much  less  likely  to 
damage  objects  [1]. 
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Figure  1 .  Electromagnetic  spectrum  showing  the  location  of  the  terahertz  band. 

From  [2] 


Therefore,  two  of  the  main  interests  for  THz  imaging  is  security  screening  [3]  and 
medical  imaging  [4],  [5],  In  the  case  of  security  screening,  metallic  objects  such  as  knives 
are  easily  observed  and  explosives  are  known  to  have  unique  absorption  spectra  in  this 
region  allowing  them  to  be  identified  inside  packages  without  damaging  the  contents.  In 
the  case  of  medical  imaging,  much  of  the  interest  lies  in  the  absorption  of  THz  by  water. 
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Skin  cancers  tend  to  have  more  blood  vessels  than  surrounding  tissue  and  thus  more 
water  content.  This  can  allow  one  to  image  cancers  with  THz  [4],  [5],  perhaps  even  in 
real-time  during  surgery.  Cultural  heritage  conservation  [6]  is  another  fascinating 
application,  where  THz  imaging  can  be  used  to  identify  art  obscured  by  other  materials 
(such  as  multiple  layers  of  paintings)  or  help  validate  the  authenticity  of  priceless  objects 
without  damaging  them.  Additionally,  quality  control  of  food  [7]  can  make  use  of  the 
same  principles  as  medical  imaging.  Food  can  be  easily  imaged  inside  packages  to  track 
how  well  it  has  dried,  or  even  in  real  time  to  track  the  efficiency  of  the  process  [7]. 
Unfortunately,  detection  in  the  THz  range  is  quite  challenging  due  to  the  low  photon 
energy  necessitating  excessive  cooling  for  optical  schemes  and  the  high  frequency 
significantly  reducing  the  effectiveness  of  electronic  schemes;  these  issues  create  the 
well-known  “Terahertz  gap”  [8].  Sensors  are  available  in  this  area,  but  are  prohibitively 
expensive  and  bulky,  largely  restricting  them  to  the  research  or  remote  sensing  domains. 

B.  TERAHERTZ  IMAGING  USING  UNCOOLED  THERMAL  DETECTORS 

Uncooled  thermal  detectors  are  a  well-established,  compact  and  inexpensive 
technology  for  imaging  in  the  infrared,  especially  in  the  form  of  microbolometer  focal 
plane  arrays  (FPAs).  Indeed,  these  sensors  naturally  evolved  in  snakes  in  the  form  of  pit 
organs  leading  to  a  high  sensitivity,  albeit  low  resolution,  infrared  imaging  system  for  8- 
12  pm  radiation  and  to  a  lesser  extent  the  3-5  pm  range  [9].  As  this  sensor  is  forced  to 
operate  at  the  temperature  of  the  snake  (which  is  the  same  as  the  thennal  background 
since  they  are  cold-blooded),  it  is  not  surprising  that  a  similar  uncooled  sensor  would 
perform  well  when  engineered  by  humans  facing  similar  limitations. 

All  thennal  sensors,  including  those  in  snakes,  work  through  a  thennally  isolated 
sensitive  area  heating  up  when  their  thermodynamic  equilibrium  is  disrupted  by  a  change 
in  incident  electromagnetic  radiation.  Thennal  sensors  differ  in  the  way  they  measure  this 
temperature  change;  microbolometers,  one  of  the  more  dominant  imaging  technologies 
for  uncooled  8-12  pm  imaging,  work  using  a  layer  of  material  in  the  thermally  isolated 
region  (microbridge)  that  experiences  changes  in  electrical  resistance  with  temperature. 
While  all  materials  change  resistance  in  response  to  a  temperature  shift,  for  certain 
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materials  such  as  titanium  (Ti)  and  Vanadium  Oxides  (VaOx),  the  magnitude  of  this 
resistance  change  is  much  greater  making  them  especially  good  choices  [10]. 

As  demonstrated  in  Figure  2,  there  is  a  large  amount  of  power  in  the  8-12  pm 
infrared  range  at  common  temperatures,  but  power  in  the  THz  region  is  comparatively 
low.  Therefore,  active  sources  are  commonly  used  as  illuminators  [3]— [7],  [1 1]— [13]  for 
THz  imaging  systems.  Possible  active  sources  include  free  electron  lasers  (FELs)  [14], 
quantum  cascade  lasers  (QCLs)  [15]— [17],  diode  based  oscillators  [18],  and  nano¬ 
electrodes  in  a  photoconductive  photomixer  [19].  QCLs,  in  particular,  are  both  compact 
and  powerful,  providing  a  useful  narrowband  source  for  imaging.  One  drawback  is  that 
QCLs  operate  at  cryogenic  temperatures. 

Wavelength  (^m) 


Figure  2.  Blackbody  distributions  for  various  temperatures  with  the  THz  band  shaded  in 
blue  and  IR  in  red.  THz  makes  up  only  a  small  fraction  of  blackbody 
emission  at  common  temperatures.  From  [20]. 

THz  images  have  been  successfully  produced  using  a  microbolometer  camera  and 
QCL  source  in  a  laboratory  setting  [1 1]— [13],  allowing  one  to  image  concealed  objects, 
as  shown  in  Figure  3(a).  However,  the  sensitivity  of  these  microbolometer  cameras  is 
limited  because  the  silicon  nitride  (SiNx)  microbridge  is  not  optimized  for  high  THz 
absorption  [21],  [22]  as  illustrated  by  the  absorption  curve  in  Figure  3(b).  Indeed, 
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uncooled  IR  microbolometers  generally  have  a  noise  equivalent  power  (NEP)  in  the  THz 
range  of  approximately  300  pW/Hz  [21]  compared  to  an  NEP  of  14  pW/Hz/2  for  the 

IR  range  [10]  implying  low  responsivity  at  THz  frequencies.  Therefore,  improving  the 
absorptivity  of  these  microbolometer  pixels  should  improve  their  performance 
considerably.  Metamaterials,  with  their  subwavelength  features,  can  possess  absorption 
characteristics  not  found  in  their  constituent  materials,  which  allows  fabrication  of 
structures  with  optimized  absorption  at  THz  frequencies.  These  structures  can  be 
incorporated  in  thermal  sensors  for  enhancing  the  sensitivity  at  THz  frequencies. 


Figure  3.  Knife  covered  in  opaque  white  plastic  tape  (a)  is  visible  in  THz. 

Absorption  (b)  of  a  1  pm  thick  Silicon  Nitride  membrane,  showing  poor 
absorption  in  THz.  (b)  From  [20]. 


C.  OVERVIEW  OF  PRESENT  RESEARCH 

In  the  following  chapters,  I  present  the  research  progression  from  designing, 
fabricating,  and  measuring  metamaterial  absorbers  alone,  integrated  into  bi-material 
sensors,  and  with  a  bolometric  layer  added  [23]-[26].  Although  each  chapter  is  largely 
based  off  a  single  published  work,  some  of  the  material  has  been  reorganized  for  a  more 
logical  presentation.  Each  phase  of  the  research  required  Finite  Element  (FE)  modeling, 
clean  room  fabrication  (performed  by  my  colleague  Dr.  Fabio  Alves),  and  finally 
measurements  in  the  laboratory  using  a  Quantum  Cascade  Laser  (QCL)  source  and  a 
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Fourier  Transform  Infrared  Spectrometer  (FTIR).  The  results  of  this  research  demonstrate 
that  metamaterials  are  a  promising  choice  for  improving  the  THz  absorption  of  thermal 
sensors  without  compromising  their  real-time  imaging  perfonnance.  They  also  offer 
guidance  on  how  to  properly  integrate  them  into  microbolometer  pixels,  allowing  one  to 
construct  highly  sensitive  THz  microbolometer  focal  plane  arrays  in  the  future. 

In  Chapter  II,  metamaterials  constructed  from  fabrication-friendly  materials  are 
investigated  and  optimized.  FE  analysis  is  used  to  guide  the  fabrication  of  metamaterial 
absorbers  optimized  to  absorb  nearly  100%  of  incoming  THz  power  at  a  specific  resonant 
frequency.  The  metamaterials  are  constructed  from  A1  and  Si-rich  SiOx;  both  of  these 
materials  are  very  suitable  for  construction  of  sensors,  especially  bi-material  sensors  as 
elaborated  on  in  Chapter  III.  The  specular  reflection  of  THz  radiation  by  these 
metamaterials  is  measured  using  our  FTIR  and  used  to  determine  the  THz  absorption  of 
these  layers  (transmission  and  higher  order  scattering  are  negligible).  These  results 
provide  us  with  the  infonnation  needed  to  tune  metamaterial  absorbing  layers  to  the 
desired  frequency  with  nearly  100%  absorption.  Finally,  the  emission  properties  are  also 
investigated  showing  that  these  materials  obey  Kirchoffs  law.  The  majority  of  the 
content  of  this  chapter  was  published  in  largely  the  same  fonn  in  Optical  Engineering 
[23],  and  presented  at  the  international  society  for  optics  and  photonics  (SPIE)  Defense 
and  Security  Screening  conference  in  April  2012.  The  emission  section  was  published  in 
essentially  the  same  form  in  Optics  Express  [24],  and  presented  at  the  SPIE  Photonics 
West  conference  in  February  2013. 

Chapter  III  uses  the  knowledge  gained  in  Chapter  II  to  design  and  integrate 
metamaterial  absorbers  into  bi-material  pixels.  First,  the  basic  physical  limitations  of 
thennal  sensors  are  explored  and  the  need  for  layers  with  a  small  thermal  capacitance  is 
established.  This  section  is  extended  from  a  paper  published  in  Optical  Materials  Express 
[25].  After  this,  the  details  behind  design,  fabrication,  and  measurement  of  bi-material 
pixels  integrated  with  THz  metamaterial  absorbers  is  explored.  FE  modeling  is  used  to 
optimize  the  thermomechanical  sensitivity  of  the  pixels,  while  the  data  from  Chapter  II  is 
used  to  optimize  absorption.  The  sensors  are  fabricated  using  standard  cleanroom 
microfabrication  techniques  and  their  performance  is  characterized  using  a  QCL  source 
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combined  with  an  optical  readout  setup  to  detennine  the  deflection  under  THz 
absorption.  This  section  demonstrates  that  the  metamaterials  in  Chapter  II  can  be 
integrated  into  bi-material  pixels,  increasing  their  THz  absorptivity.  Other  than  the  first 
section,  this  work  is  contained  entirely  within  a  paper  published  in  Optics  Express  [26], 
and  was  also  presented  at  the  SPIE  Photonics  West  conference  in  February  2013. 

Although  bi-material  pixels  face  the  same  thermal  challenges  that  microbolometer 
pixels  do,  suggesting  that  metamaterials  will  enhance  their  performance,  they  do  not  have 
the  same  structure.  Indeed,  the  effect  of  the  bolometric  layer  in  a  microbolometer  pixel 
must  be  explored  to  ensure  that  metamaterials  can  incorporate  this  layer  without  having 
their  absorption  adversely  affected,  or  the  thermal  properties  of  the  pixels  compromised. 
Therefore,  to  investigate  this  metamaterials  with  an  embedded  bolometric  layer  (in  this 
case,  Ti)  were  fabricated.  These  metamaterials  were  fabricated  with  the  bolometric  layer 
at  several  different  depths  and  the  subsequent  FTIR  measurements  demonstrated  that  if 
the  bolometric  layer  is  close  to  the  ground  plane,  the  absorption  is  largely  unaffected. 
This  information  was  used  to  validate  the  FE  model,  which  also  suggests  that  a 
bolometric  layer  with  a  low  sheet  resistance  will  also  have  a  smaller  impact  on  the 
maximum  absorption  amplitude.  Finally,  different  possible  configurations  are  explored 
for  their  various  advantages  and  disadvantages  and  future  work  in  the  area  is  suggested. 
Other  than  the  last  sections,  this  chapter  is  published  practically  unaltered  in  Optical 
Materials  Express  [25]. 

The  progression  from  metamaterial  absorber  design,  integration  into  thermal 
sensors,  and  analysis  of  the  impact  of  a  bolometric  layer  represent  significant  progress 
towards  the  successful  integration  of  THz  metamaterial  absorbers  to  enhance  the 
performance  of  thennal  sensors,  specifically  microbolometers.  Chapter  II  provides  the 
general  information  necessary  when  constructing  these  layers.  The  results  from  Chapter 
III  confirm  that  these  metamaterials  can  be  integrated  into  thennal  sensors  and  that 
constituent  materials  that  maximize  the  thermal  sensitivity  of  the  sensor  can  be  used 
while  still  maintaining  high  THz  absorptivity  and  a  small  thickness,  highlighting  the 
primary  advantage  of  metamaterials.  Finally,  Chapter  IV  confirms  that  a  bolometric  layer 
can  be  integrated  and  suggests  that  this  layer  should  be  placed  as  close  as  possible  to  the 
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ground  plane  and  have  as  small  of  a  sheet  resistance  as  possible.  Ideally,  these  advances 
will  lead  to  highly  sensitive  THz  microbolometer  cameras  using  metamaterial  absorbers 
for  practical  real-time  imaging  applications. 
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II.  METAMATERIAL  ABSORBERS 


Note:  The  majority  of  the  content  of  this  chapter  was  published  in  largely  the 
same  form  in  Optical  Engineering  [23],  and  presented  at  the  international  society  for 
optics  and  photonics  (SPIE)  Defense  and  Security  Screening  Conference  in  April  2012. 
The  emission  section  was  published  in  essentially  the  same  form  in  Optics  Express  [24], 
and  presented  at  the  SPIE  Photonics  West  Conference  in  February  2013. 

Metamaterials  are  artificial  structures  with  subwavelength  features  that  can 
exhibit  overall  optical  properties  not  characteristic  of  their  constituent  materials.  Indeed, 
they  can  even  have  characteristics  that  are  not  found  in  nature,  such  as  a  negative 
refractive  index  [27].  Absorption  in  these  structures,  especially  through  resistive  losses, 
can  be  a  vexing  problem  for  applications  such  as  superlenses  [28].  However,  certain 
applications  can  benefit  greatly  from  enhanced  absorption,  such  as  thennal  sensors. 
Landy  et  al.  [29]  initially  envisioned  a  material  which  would  both  not  transmit  light  yet 
would  be  impedance  matched  to  free  space  through  adjustment  of  the  dielectric  and 
magnetic  response  of  the  metamaterial.  This  would  result  in  a  material  that  could  not 
reflect  or  transmit  light,  and  therefore  must  absorb  all  incident  radiation.  In  practice,  this 
typically  is  only  achievable  over  certain  frequency  bands,  depending  on  the  geometry  of 
the  subwavelength  features.  Materials  with  this  property  have  been  realized  across  the 
electromagnetic  spectrum  [30]— [43].  In  principle,  one  can  fabricate  a  metamaterial  for 
any  frequency  within  this  range  of  interest  by  an  appropriate  adjustment  to  the  geometry 
of  the  subwavelength  features.  Fortunately,  this  tuning  is  only  weakly  dependent  on  the 
materials  used,  so  these  metamaterial  absorbing  films  can  be  constructed  using  standard, 
mature  MEMS  microfabrication  techniques.  These  metamaterial  films  typically  consist  of 
a  periodic-patterned  metal  layer,  a  dielectric  spacer,  and  a  ground  plane  [30]— [43].  In  this 
work,  a  periodic  array  of  metal  squares  was  used  as  shown  in  Figure  4.  Metal  squares 
were  chosen  due  to  the  ease  of  their  fabrication  and  thin  dielectric  spacer  requirements. 
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A. 


SIMULATION 


While  empirical  studies  are,  of  course,  essential  to  evaluate  the  properties  of 
metamaterials,  numerical  modeling  can  very  effectively  predict  the  properties  of  a 
metamaterial  absorber  without  the  need  to  fabricate  every  possible  geometry.  Reflection 
is  suppressed  as  incident  EM  waves  excite  localized  surface  plasmon  resonances,  leading 
to  resistive  losses  in  the  metal  and  additional  losses  in  the  dielectric  layers  converting 
some  of  the  incident  energy  to  heat,  allowing  the  incident  radiation  to  be  absorbed  [30], 
[32].  Resonant  absorption  frequency  (v)  and  absorption  magnitude  can  be  understood 
using  explanations  such  as  impedance-matching  to  free  space  [29],  an  equivalent  resistor 
inductor  capacitor  (RLC)  circuit  [30],  [32],  confined  TM  cavity  modes  [44],  interference 
of  multiple  reflections  [45],  and  transmission  lines  [46].  However,  due  to  the  complex 
nature  of  metamaterials,  these  models  are  typically  guided  by  finite  element  (FE) 
modeling  (of  at  least  the  patterned  metallic  layer)  [30],  [32],  [44]-[46],  Therefore, 
understanding  numerical  models  of  metamaterial  behavior  is  essential  to  predicting  their 
properties.  In  this  work,  I  used  COMSOL  Finite  Element  modeling  software.  COMSOL 
is  a  multiphysics  package  with  many  different  modules  available,  but  the  one  used  for 
analyzing  metamaterials  is  the  radio  frequency  (RF)  module.  Great  care  must  be  taken 
when  constructing  3D  models,  so  that  one  can  preserve  the  physics  while  still  obtaining  a 
result  in  a  reasonable  amount  of  time. 

When  creating  a  FE  model,  boundary  conditions  and  meshing  are  crucial 
considerations.  The  geometry  itself  is  usually  fairly  straightforward  to  create  and  so  is  of 
little  concern.  I  should  emphasize  that  the  geometry  should  match  the  fabricated  structure 
as  accurately  as  possible,  as  even  minor  changes  in  geometry  can  lead  to  variations  in 
absorption  characteristics.  A  typical  rule  of  thumb  for  FE  analysis  is  that  elements  should 
be  no  further  than  A/5  apart,  which  can  result  in  an  excessive  amount  of  mesh  if  a  very 
large  volume  of  metamaterial  is  simulated  in  a  3D  model.  To  simplify  the  situation,  we 
can  exploit  the  fact  that  metamaterials  are  periodic  structures  and  thus  any 
electromagnetic  solution  must  satisfy  periodic  boundary  conditions.  Therefore,  only  a 
single  unit  cell  of  metamaterial  (see  Figure  4(a)),  with  appropriate  periodic  boundary 
conditions,  is  necessary  for  FE  analysis.  For  radiation  normally  incident  to  the  surface, 

10 


the  electric  and  magnetic  fields  of  an  incident  plane  wave  can  be  broken  into  two 
polarizations  with  fields  normal  to  the  boundaries.  In  this  case,  a  perfect  electric 
conductor  (PEC)  will  enforce  periodic  boundary  conditions  for  electric  fields  normal  to 
the  boundary,  while  a  perfect  magnetic  conductor  (PMC)  will  do  the  same  for  magnetic 
fields  normal  to  a  boundary.  In  the  case  of  radiation  not  normally  incident  to  a  surface, 
Floquet  boundary  conditions  are  appropriate  as  they  attach  the  necessary  phase  factor  to 
the  radiation  [47].  However,  in  most  thermal  sensor  applications  radiation  is  near  normal 
incidence  making  normally  incident  radiation  a  good  approximation. 


Scattering  Boundary  Condition 

(a)  (b)  (c) 

Figure  4.  Unit  cell  (a)  of  a  metamaterial  with  a  single  absorption  peak.  The  model  also 
incorporates  the  unintentional  over-etch  into  the  SiOx  during  the  fabrication. 
Schematic  representation  (b)  of  the  unit  cell  in  the  FE  modeling  program  with 
boundary  conditions  and  THz  plane  wave  at  normal  incidence.  Time- 
averaged  power  flow  (c)  near  the  resonant  absorption  frequency  and  arrow 
lengths  relate  to  magnitude  in  log  scale.  From  [23] 


The  model  must  also  allow  electromagnetic  radiation  to  interact  with  the 
metamaterial  absorbers.  To  accomplish  this,  one  must  choose  the  proper  boundary 
conditions  for  the  surfaces  intended  to  allow  radiation  to  penetrate  through.  One  surface 
must  somehow  generate,  or  allow  penetration  of,  incident  radiation  and  also  allow 
reflected  radiation  to  leave  the  model  or  absorb  it  without  allowing  any  radiation  to 
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reflect  back  in.  The  second  surface,  typically  “below”  the  ground  plane,  must  allow  any 
transmitted  radiation  to  leave  the  model  or  absorb  it,  also  without  any  radiation  reflecting 
back  into  the  model.  Several  popular  methods  available  in  COMSOL  are  scattering 
boundary  conditions,  internal  ports,  and  perfectly  matched  layers  (PMLs).  Scattering 
boundary  conditions  are  meant  to  approximate  an  infinite  space  or  “open”  boundary  for 
radiation  to  propagate  out  into.  They  work  as  first  order  absorbing  boundary  conditions 
for  plane,  cylindrical,  or  spherical  waves.  Port  boundary  conditions  are  similar  in  that 
they  are  perfectly  absorbing  boundary  conditions  for  general  modes  of  known  shape. 
Both  of  these  boundaries  can  also  allow  the  same  types  of  waves  to  penetrate  the 
boundary  into  the  model,  supplying  the  incident  radiation.  The  power  is  easily  controlled 

by  adjusting  the  Poynting  vector  and  can  be  calculated  using  A^S-n^J  for  plane 

boundaries.  While  not  exactly  a  boundary  condition,  a  PML  can  also  be  used  to 
approximate  an  open  boundary.  A  PML  is  actually  very  similar  to  the  initial  idea  of  a 
metamaterial  absorber;  it  is  impedance  matched  to  the  material  that  it  is  adjacent  to  so 
that  it  can  eliminate  reflection  and  uses  a  complex  coordinate  transformation  to  cause  an 
incident  electromagnetic  wave  to  decay  as  it  propagates  resulting  in  essentially  total 
absorption  [47].  The  PML  actually  has  a  finite  thickness  and  must  have  additional  mesh 
elements,  and  should  usually  be  coupled  with  a  scattering  boundary  condition  to 
minimize  back  reflections.  All  of  these  tools  tend  to  work  very  well  for  a  particular  wave 
mode  and  propagation  vector  and  then  start  to  allow  some  reflection  back  into  the  model, 
so  care  should  be  taken  in  their  implementation.  A  schematic  of  this  setup  is  shown  in 
Figure  4(b). 


Once  boundary  conditions,  geometry,  and  meshing  are  finished,  the  wave 
equations  for  the  model  to  solve  must  be  specified.  The  default  option  for  COMSOL  is 
“relative  permittivity”  which  solves  the  EM  wave  equation  of  the  form: 


V  x  ju~l  (Vxl)  - kl 


- 


cos, 


E  =  0 


0  J 


(1) 


Thus,  this  equation  admits  relative  permeability  and  permittivity  as  well  as 


conductivity.  This  is  convenient  for  materials  that  can  be  defined  entirely  by  permeability 
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and  pennittivity,  such  as  free  space,  or  those  where  conductivity  dominates,  such  as  most 
metals  in  THz.  For  some  metals,  frequency-dependent  conductivity  is  needed  if  it  varies 
appreciably  within  the  frequency  range  of  interest  (such  as  Au  [44]).  In  the  case  of  A1 
used  in  our  metamaterial  structures,  the  relaxation  is  fast  enough  that  the  Drude  model 
approximates  to  the  DC  conductivity  [49].  It  should  be  noted  that  thin  films  of  metals  do 
not  have  the  same  conductivity  as  bulk  metals  and  also  may  have  variable  conductivities 
depending  on  deposition  technique  and  surface  [49].  COMSOL  does  have  a  “Drude- 
Lorentz  dispersion  model”  option,  but  this  approximation  renders  it  unnecessary  for  A1  as 
absorption  remains  largely  unchanged.  The  second  convenient  setting  is  “Refractive 
index”  which  solves  the  EM  wave  equation  of  the  form: 

V  x  (V  xfi'j-  k\srE  =  0  (2) 

Clearly,  this  form  assumes  relative  permeability  is  1  and  conductivity  is  contained 
in  the  relative  pennittivity  via  the  extinction  coefficient  [46],  It  then  calculates  the 

pennittivity  using  the  complex  refractive  index:  er=(n-ik\ ,  using  the  usual 

convention  that  a  positive  “k”  refers  to  loss  while  a  negative  “k”  is  gain.  This  fonn  is 
useful  for  dielectrics,  where  the  permeability  can  largely  be  ignored,  and  the  literature 
commonly  specifies  the  complex  refractive  index. 

Once  the  model  construction  is  complete,  the  solvers  must  be  chosen.  The  models 
in  this  work  use  a  frequency  domain  solver,  which  simply  solves  the  EM  wave  equations 
at  specified  frequencies.  This  results  in  the  E  and  B-field  solutions  allowing  one  to 
detennine  information  such  as  power  flow  (see  Figure  4(c)).  More  resolution  simply 
correlates  to  a  longer  computation  time,  but  more  precise  absorption  curves.  Of  greater 
consequence  is  the  type  of  solver  used.  Direct  solvers  are  preferable  if  the  modeling 
computer  possesses  sufficient  memory  as  there  are  no  convergence  issues,  yet  indirect 
solvers  can  solve  denser  meshes  with  the  same  resources  but  can  have  convergence 
issues.  Thus,  the  choice  of  solver  and  settings  will  largely  be  determined  on  a  case-by- 
case  basis.  Direct  solvers  were  used  for  the  models  in  this  work,  due  to  the  computational 
resources  available.  In  addition,  parametric  sweeps  can  be  used  to  run  several  different 
configurations  of  metamaterials  together. 
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Finally,  once  the  model  has  finished  its  run,  one  must  extract  absorption 
characteristics.  In  all  cases,  when  magnetic  losses  are  largely  ignorable  in  the  constituent 
materials,  a  convenient  way  to  extract  absorption  is  through  performing  volume 
integrations  of  the  resistive  losses  in  the  metamaterial.  Non-linear  effects  are  largely 
ignored  in  these  models,  so  the  incoming  power  can  be  set  to  1  W  and  the  absorption  can 
simply  be  read  off  as  the  Watts  of  power  absorbed  in  the  metamaterial.  This  also  has  the 
advantage  of  identifying  which  layers  are  actually  heating  up,  allowing  one  to  more 
effectively  model  metamaterials  that  are  not  thermally  linked  as  in  Figure  5. 


Figure  5.  Modeled  absorption  characteristics  of  the  three  layers  of  a  metamaterial  with 

single  absorption  peak.  From  [23] 


Additionally,  when  using  ports,  one  can  also  extract  transmission  and  reflection 
coefficients  by  inverting  the  S  parameters.  This  method  is  also  important  for  the 
reflection  model  that  relies  on  the  resonator  layer  having  some  frequency  dependent 
reflection  coefficient  [45].  This,  of  course,  assumes  that  higher  order  scattering  is  largely 
ignorable  which  should  be  true  for  sub-wavelength  features  of  this  kind  and  is  known  to 
be  small  for  metamaterials  in  other  frequency  ranges  [46].  In  practice,  transmission 
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through  a  metamaterial  is  vanishingly  small  due  to  the  optically  thick  ground  plane  of 
approximately  100  mn.  The  skin  depth  of  A1  layers  is  approximately  80  mn  for  1  THz, 
and  both  models  and  measurements  show  no  significant  transmission.  Finally,  reflection 
and  transmission  coefficient  techniques  and  resistive  loss  can  be  used  together  to  validate 
that  models  are  working  correctly  [30], 

B.  FABRICATION 

Metamaterial  absorbers,  whether  split  ring  resonators  or  arrays  of  squares, 
typically  consist  of  a  patterned  metal  layer,  a  dielectric  spacer  of  some  kind,  and  a  ground 
plane.  Variations  including  a  bolometric  layer  and  integration  into  bi-material  sensors 
will  be  discussed  in  later  chapters.  Fabrication  of  these  layers  is  easily  done  using 
standard  microfabrication  techniques.  This  was  performed  at  Oak  Ridge  National 
Laboratory  by  Doctor  Fabio  Alves  and  Assistant  Professor  Dragoslav  Grbovic.  At  every 
point  in  the  fabrication,  the  materials  were  carefully  monitored  for  relevant  parameters 
such  as  residual  stress  and  DC  conductivity.  These  measurements  were  perfonned  with 
stylus  profilometry  and  optical  interferometry  in  the  case  of  thickness,  four-point  probe 
measurements  in  the  case  of  conductivity,  and  wafer  curvature  measurements  using  a 
stylus  or  optical  profilometer  to  evaluate  residual  stress.  As  transmission  is  essentially 
zero  due  to  the  thick  ground  plane  of  the  metamaterial  [30],  [48],  one  does  not  need  to 
use  low  conductivity  substrates  as  with  thin  metal  films  [49].  If  conductivity  of  the  wafer 
is  a  concern,  a  layer  of  insulating  dielectric  can  be  deposited.  Therefore,  a  single  side 
polished  (SSP)  Si  wafer  of  any  thickness  or  doping  concentration  should  be  suitable,  but 
thicker  wafers  are  less  fragile.  Generally  500  pm  thick  wafers  were  used  due  to  the  low 
price  and  durability  of  these  wafers. 

The  first  step  is  to  deposit  the  metallic  ground  plane.  We  typically  used  A1 
deposited  using  an  e-beam  evaporator  due  to  the  fast  deposition  rate,  good  adhesion,  and 
good  conductivity.  Typically,  these  A1  layers  were  measured  to  have  a  DC  conductivity 
of  107  Siemens/meter  for  layers  of  approximately  100  mn.  The  bulk  conductivity  of  A1  is 
over  3  times  larger  than  that  of  the  A1  films,  therefore  measurement  with  a  four-point 
probe  is  essential.  The  second  step  is  the  dielectric  spacer,  which  in  this  work  was  SiOx. 
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Specifically,  this  was  a  low  stress  oxide  deposited  with  a  Plasma  Enhanced  Chemical 
Vapor  Deposition  (PECVD)  system  at  350°C  with  more  Si  than  the  standard  SiC>2 
formulation  (-13  MPa  compared  to  -140  MPa).  Dielectric  film  parameters  can  vary 
between  machines  and  even  through  small  changes  in  machine  parameters;  therefore, 
great  care  was  taken  to  reproduce  the  same  SiOx  for  every  phase  of  the  project.  Finally,  a 
second  A1  layer  is  deposited  on  top  of  this  dielectric.  The  layer  is  deposited  using  exactly 
the  same  settings  as  the  ground  plane  metal.  To  pattern  this  layer  and  form  a 
metamaterial,  it  must  be  etched.  First,  after  the  wafer  is  prepared,  photoresist  is  spun  on 
using  a  spin  coater.  After  some  baking  steps  specific  to  the  photoresist,  the  wafer  is 
placed  into  a  contact  aligner.  A  mask  with  the  metamaterial  absorber  resonators  is  then 
placed  on  top  of  the  wafer  and  the  photoresist  is  exposed  to  UV  light.  After  this,  the 
photoresist  is  developed  and  processed  further,  again  with  steps  specific  to  the  type  of 
photoresist.  In  this  work,  positive  photoresist  was  used,  so  unexposed  areas  remained  on 
the  wafer  while  exposed  areas  were  dissolved.  To  remove  the  unwanted,  now  exposed  Al, 
an  Ar  plasma  sputter  etcher  was  used.  The  sputter  etcher  accelerates  Ar  ions  at  the  wafer, 
removing  both  photoresist  and  exposed  Al.  This  etcher  is  anisotropic  (directional)  which 
reduces  the  “undercut”  one  would  get  from  isotropic  etching  techniques  such  as  acid 
etching.  The  remaining  photoresist  protects  the  Al  that  we  wish  to  stay  on  the  wafer. 
Unfortunately,  this  process  is  not  very  selective  between  Al  and  the  underlying  dielectric 
and  the  etching  must  occur  long  enough  to  go  all  the  way  through  the  Al  layer.  This 
results  in  an  effect  known  as  “over-etch”  where  some  of  the  dielectric  is  also  removed  as 
illustrated  in  Figure  4(a).  The  effect  of  over-etch  can  be  accounted  for  in  simulations  by 
simply  adding  to  the  geometry.  Finally,  the  photoresist  is  removed  using  appropriate 
chemicals  and  an  isotropic  oxygen  plasma  “barrel”  etcher.  The  result  is  a  wafer  with  25 
different  metamaterials  pictured  in  Figure  6. 

It  should  be  noted  that  other  materials  and  processes  can  be  used.  Indeed,  acid 
etching  (aka  “wet”  etching)  is  much  more  selective  and  could  result  in  no  over-etch  into 
the  dielectric.  If  the  etching  is  carefully  controlled  and  reproducible,  it  may  well  be  a 
better  choice  than  sputter  etching.  Other  suitable  choices  for  the  Al  etching  include 
reactive  ion  etching  and  liftoff.  The  advantage  of  reactive  ion  etching  is  its  better 
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selectivity  towards  Al.  In  the  case  of  liftoff,  the  photoresist  is  deposited  first.  Inverted 
masks  are  used  such  that  the  unwanted  Al  is  deposited  on  the  photoresist  and  is  removed 
with  it.  This  requires  no  metal  etchant  at  all,  and  is  therefore  faster  and  easier,  albeit  with 
less  precision.  Also,  other  materials  can  be  used  for  metals  include  Au  or  Ag  or 
dielectrics  such  as  SiNx  and  polyimides. 


Figure  6.  One  of  the  fabricated  wafers  with  25  different  metamaterials  arranged  radially 
for  automated  reflection  measurements  using  a  PIKE  MappIR  FTIR 

accessory.  From  [23] 


C.  MEASUREMENT 

To  evaluate  the  THz  optical  properties  of  the  metamaterial  structures,  a  Fourier 
Transform  Infrared  Spectrometer  (FTIR)  was  used.  This  device  consists  of  a  Michelson 
interferometer  with  a  moving  mirror  in  one  arm;  this  results  in  periodic  oscillations  in 
amplitude  whose  frequency  depends  on  the  wavelength  of  radiation  and  velocity  of  the 
mirror.  Thus,  the  FTIR’s  thermal  detector  can  measure  every  frequency  simultaneously 
and  extract  the  frequency  data  with  a  Discrete  Fourier  Transform  (DFT).  This  powerful 
technique  is  applicable  across  the  frequency  spectrum,  but  is  particularly  effective  for  the 
Infrared  region  of  the  spectrum.  The  FTIR  available  at  NPS  is  a  Thermo-Nicolet  Nexus 
870  fitted  with  a  globar  source  (hot  SiC  filament)  and  a  Deuterated  Triglycine  Sulfate 
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(DTGS)  pyroelectric  detector  with  a  polyethelene  window.  This  is  the  recommended 
configuration  for  Far-IR  spectrometry  (defined  by  Nicolet  as  50-600  cm'1  or  1.5-30 
THz);  however,  the  1.5-3  THz  range  is  difficult  to  use  due  to  the  fact  that  the  source 
power  is  relatively  low  as  shown  in  Figure  7.  For  measurements  in  the  Mid  Infrared 
(Mid-IR)  range  (defined  by  Nicolet  as  7400-350  cm’1  or  1.35-28.6  pm),  the  beam  splitter 
is  changed  to  KBr  and  the  detector  is  changed  to  a  Thennoelectrically  cooled  DTGS 
detector  with  a  KBr  window.  For  metamaterial  absorption,  the  Far-IR  configuration  is  of 
important  as  it  contains  the  THz  range  of  interest.  In  terms  of  FTIR  settings  and 
configurations,  one  must  achieve  a  large  enough  signal  to  align  the  FTIR,  yet  a  small 
enough  signal  to  avoid  saturation  of  the  detector.  In  the  default  configuration  of  the  FTIR, 
it  can  only  measure  transmission  through  a  sample. 


Figure  7.  Raw  FTIR  data  for  a  metamaterial  absorber,  showing  the  excessive  noise  in 

the  spectral  region  below  3  THz. 


Clever  use  of  optics  can  allow  one  to  measure  reflection  as  shown  in  Figure  8; 
specifically,  an  FTIR  accessory  known  as  MappIR  can  be  programmed  to  automatically 
measure  either  transmission  or  reflection  at  predefined  points  on  a  wafer  resulting  in  a 
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significant  time  savings  when  dealing  with  a  large  number  of  metamaterials.  In  the  case 
of  emission  measurements,  the  FTIR  is  programmed  to  see  the  metamaterials  and  a  black 
carbon  blackbody  (prepared  by  depositing  candle  soot  on  a  suitable  substrate)  as  an 
external  source.  An  off  axis  parabolic  mirror  collimates  the  thermal  radiation  from  a 
small  area  of  the  sample  on  a  hot  plate  and  directs  it  into  the  FTIR.  In  all  of  these 
configurations,  the  THz  must  pass  through  the  outside  air,  leading  to  frequency 
dependent  losses  that  cannot  be  ignored.  The  inner  area  of  the  FTIR  is  purged  with  clean, 
dry,  low  CCF  air  to  preserve  the  optics  and  boost  signal  strength.  Despite  this  a 
background  spectrum  must  always  be  taken  in  FTIR  measurements  to  subtract  the  non- 
uniform  spectral  characteristics  of  the  system.  Appropriate  backgrounds  depend  on  the 
measurement  and  configuration,  but  typically  are  air  for  transmission,  thick  uniform  Au 
or  A1  films  for  reflection,  or  black  carbon  blackbodies  for  emission.  The  FTIR  can 
automatically  compare  these  measurements  to  the  sample  or  they  can  be  done  in  the 
graphing  software.  Typically  the  output  is  given  in  “percent  transmission”  which 
essentially  divides  the  sample  points  by  the  background  points.  This  results  in  sizable 
errors  when  the  signal  is  much  smaller  than  the  noise  as  seen  for  the  1.5-3  TFIz  range  in 
Figure  7. 


(a)  (b) 


Figure  8.  Schematic  representation  of  the  experimental  configuration  used  for 

measuring  the  reflection  (a)  and  transmission  (b),  before  using  the  MappIR 
accessory.  The  accessory  performs  essentially  the  same  measurements,  but 
operates  at  an  incident  angle  of  15  degrees.  From  [49]. 
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D. 


OPTIMIZATION  AND  TUNING  OF  SINGLE-BAND  METAMATERIAL 
ABSORBERS 


The  simplest  metamaterial  geometry  contains  only  one  type  of  square  arranged  in 
a  periodic  array.  I  refer  to  them  as  “single-band”  absorbers  because  they  have  only  one 
absorption  peak  in  the  frequency  range  of  interest.  As  Figure  9  demonstrates,  we  see  that 
the  lowest  order  confined  mode  consists  of  a  single  oscillating  dipole  formed  by  the 
incident  radiation  interacting  with  the  square  with  opposing  currents  and  charges  in  the 
ground  plane. 


Figure  9.  FE  model  solution  for  16  pm  squares  at  the  resonant  absorption  frequency. 

Electric  (red)  and  magnetic  (blue)  fields  for  a  metamaterial  absorber  at  its 
resonant  frequency.  Note  the  opposing  electric  fields  on  either  side  of  the 
square  and  the  reversal  of  the  magnetic  field. 


This  suggests  that  we  should  expect  the  wavelength  of  any  trapped  first-order 
modes  should  be  proportional  to  the  size  of  the  structure  in  the  direction  of  the  incident 
E-field  component  parallel  to  the  surface.  Of  course,  this  does  not  take  into  account  the 
boundary  conditions  or  the  coupling  between  squares.  Figure  10(a)  shows  experimental 
results  for  single  band  absorbers  with  the  same  21  pm  periodicity,  but  incrementally 
increasing  square  size.  As  expected,  the  resonant  frequency  decreases  as  square  size 
increases.  Figure  10(b)  shows  that  the  resonant  frequency  is  roughly  proportional  to  the 
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inverse  of  the  square  size,  as  expected  due  the  fact  that  wavelength  is  inversely 
proportional  to  frequency.  Therefore,  modifying  square  size  is  an  effective  method  for 
controlling  the  resonant  absorption  frequency. 


Frequency  (THz)  Inverse  Square  Size  (nm'1) 
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Figure  10.  FTIR  measurements  (a)  for  varying  square  size,  21  pm  pitch,  and  1.47  pm 
SiOx  .  Linear  fit  (b)  of  inverse  square  size  (1/s)  to  resonant  absorption 
frequency  (v).  The  effects  of  varying  dielectric  thickness  (c)  for  12  pm 
squares  and  21  pm  pitch  and  varying  pitch  (d)  with  16  pm  squares  and  1.47 

pm  SiOx.  From  [23] 


Now  that  the  most  effective  frequency  tuning  method  has  been  established,  we  must 
optimize  the  remaining  geometric  parameters  to  achieve  the  highest  possible  maximum 
absorption  amplitude.  Dielectric  thickness  and  periodicity  are  known  to  be  significant 
parameters  for  maximum  absorption.  As  Figure  10(c)  demonstrates,  a  higher  dielectric 
thicknesses  typically  result  in  higher  absorption.  As  no  reduction  in  absorption  amplitude 
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was  observed  for  increased  thickness  in  this  range,  it  seems  that  all  absorbers  have  an 
optimal  thickness  at  or  above  1.66  pm.  Once  past  this  optimal  thickness,  absorption 
magnitude  will  decrease  with  increasing  thickness  until  another  optimal  thickness  is 
approached.  The  reflection  model  of  absorption  suggests  that,  similar  to  thin  film 
interference,  there  should  be  several  optimal  thicknesses.  In  the  case  of  thermal  sensor 
applications,  we  will  see  that  in  most  cases  these  materials  should  be  as  thin  as  possible. 
Periodicity,  or  pitch,  shows  a  similar  relationship  to  dielectric  thickness,  showing  more 
significant  frequency  shifts  as  demonstrated  in  Figure  10(d).  Intuitively,  a  smaller  pitch 
should  correspond  to  a  higher  absorption  due  to  a  higher  density  of  squares  [30],  [32]  or 
better  coupling  to  external  radiation  [44],  It  should  be  noted  that  there  is  a  relationship 
between  pitch  and  dielectric  thickness  such  that  a  given  dielectric  thickness  will  have  an 
optimal  pitch.  As  Figure  1 1  shows,  simulations  suggest  that  at  certain  thicknesses 
increasing  pitch  can  increase  absorption  magnitude.  While  this  is  certainly  relevant  if  one 
wishes  to  construct  very  thick  absorbers,  in  the  thin  regime  a  decrease  in  pitch  will 
always  be  desirable.  Of  course,  any  resonant  absorption  frequency  shifts  due  to 
periodicity  or  dielectric  thickness  changes  must  be  compensated  for  through  tuning  the 
square  size. 

While  metamaterial  absorbers  can  be  constructed  using  many  different  materials, 
the  materials  used  do  have  a  significant  impact  on  the  optimal  dimensions  for  the 
metamaterial  geometry.  Metals,  or  equivalent  high  conductivity  materials,  must  have 
sufficient  thickness  and  conductivity  to  function  as  resonators.  Additionally,  some  metals 
such  as  Au  have  conductivities  that  vary  greatly  across  1-  10  THz  [44],  [45],  [48]. 
However,  the  A1  films  used  in  this  experiment  have  essentially  the  same  conductivity 
across  the  3-8  THz  range  [48].  The  refractive  index  of  the  dielectric  spacer  is  also 
significant,  as  it  will  relate  to  the  capacitance  of  the  squares  or  the  wavelength  of  the 
trapped  dipole  modes,  resulting  in  higher  absorption  frequencies  for  lower  refractive 
indexes.  Additionally,  processing  steps  that  over-etch  into  the  dielectric  cannot  be 
ignored.  This  removes  dielectric  from  around  the  squares,  causing  a  shift  in  resonant 
absorption  frequency.  This  could,  of  course,  be  used  intentionally  to  shift  resonant 
frequency  and  reduce  the  thermal  mass  of  the  absorbers  if  desired. 
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Figure  1 1 .  Simulations  for  1 6  pm  square  size  and  2 1  pm  pitch  with  two  dielectric 

thicknesses  showing  that  for  thicker  cases,  increasing  periodicity  can,  in  fact, 

increase  percent  absorption. 


E.  OTHER  METAMATERIAL  CONFIGURATIONS  AND  EMISSION 

CHARACTERISTICS 

While  single  band  metamaterial  absorbers  are  ideal  for  coupling  with  a  single 
narrowband  source  such  as  a  QCL  for  active  THZ  imaging,  in  some  cases  a  broader 
absorption  peak  may  be  desired.  Additionally  one  might  require  two  absorption  peaks,  for 
instance  when  working  with  two  separate  QCL  sources.  Metamaterials  can  be  fabricated 
with  these  properties  by  using  more  complex  unit  cells.  A  possible  configuration  is  to 
simply  alternate  the  squares  resulting  in  a  pattern  reminiscent  of  a  checkerboard,  as 
illustrated  in  Figure  12(a)  with  squares  dimensions  of  13  and  17  pm  and  a  dielectric  layer 
thickness  of  1.47  pm.  Figure  12  shows  the  measured  (b)  and  simulated  (c)  absorption  of 
structure  in  (a)  along  with  absorption  spectra  by  removing  one  type  of  squares  and 
keeping  the  rest  of  the  parameters  the  same. 
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Figure  12.  Unit  cell  (a)  of  a  dual  band  metamaterial  configuration.  Measured  (b)  and 
modeled  (c)  absorption  for  a  dual-band  absorber  with  squares  that  have  a  4 
pm  size  difference  compared  to  the  same  configuration  with  only  one  type  of 
square.  Comparison  (d)  of  the  absorption  spectrum  for  dual  band  absorber 
with  2  pm  square  size  difference,  overall  pitch  of  42  pm  and  SiOx  thickness 
1 .47  pm  to  a  single  band  absorber  with  an  intermediate  square  size.  From 

[23]. 


It  can  be  easily  seen  from  Figure  12(b)  that  each  type  of  square  is  responsible  for 
one  of  the  two  peaks  of  the  combined  structure  absorption.  For  the  configuration  shown 
in  Figure  12(a),  the  absorption  measurement  shows  two  well-separated  absorption  peaks 
at  4.25  and  5.5  THz  similar  to  that  obtained  in  the  simulation.  In  addition,  Figure  12(b) 
shows  measured  single-band  absorption  from  two  complementary  structures  fabricated 
with  only  one  type  of  square  while  maintaining  the  same  dimension  and  pitch.  In  this 
case,  single-absorption  peaks  clearly  match  with  the  absorption  peaks  of  the  dual-band 
absorber,  indicating  that  each  set  of  squares  in  Figure  12(a)  is  responsible  for  generating 
the  two  peaks  as  predicted  by  the  FE  model  (Figure  12(c)).  As  the  difference  in 
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dimensions  of  square  size  is  reduced  (16  and  18  pm),  the  absorption  peaks  begin  to 
overlap  and  merge  together  as  illustrated  by  measured  data  in  Figure  12(d).  Thus,  by 
bringing  the  size  of  the  squares  closer  together,  it  is  possible  to  produce  a  broader 
absorption  peak  than  a  single  square  configuration  of  intermediate  size  (17  pm).  In  both 
the  measurements  and  the  FE  model  (Figure  12  (b)  and  (c),  respectively),  it  is  apparent 
that  the  lower  frequency  absorption  is  enhanced  while  the  higher  frequency  peak  is 
suppressed.  This  causes  the  asymmetric  peak  shown  in  Figure  12(d),  as  the  two 
absorption  peaks  are  merged  by  making  the  dimension  of  squares  close  to  each  other. 
This  can  partially  be  explained  by  the  fact  that  the  smaller  squares  are  more  active  at  the 
larger  square  resonant  frequency  than  vice  versa  as  shown  in  Figure  13. 


Figure  13.  FE  model  solution  for  16  and  18  pm  squares  at  the  resonant  absorption 
frequency  of  the  18  pm  squares  (a)  and  the  16  pm  squares  (b).  The  center 
square  is  18  pm.  Electric  (red)  and  magnetic  (blue)  fields  are  shown.  Note 
how  the  16  pm  squares  have  a  larger  electric  field  in  (a)  than  the  18  pm 
squares  do  in  (b)  showing  that  the  16  pm  squares  assist  the  18  pm  at  their 
resonance,  but  not  vice  versa. 


Another  configuration  for  producing  two  absorption  peaks  is  shown  in  Figure 


14(a)  where  smaller  squares  are  placed  around  larger  ones  in  a  tile-like  pattern.  This 
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configuration  allows  for  squares  with  large  size  differences,  and  therefore  absorption 
peaks  at  widely  separated  frequencies,  to  be  more  densely  packed  to  increase  absorption 
magnitude  as  demonstrated  in  Figure  14(b).  Measurements  showed  two  distinct 
absorption  peaks  at  frequencies  4  and  6.5  THz,  as  depicted  in  Figure  14(b),  with  the  high 
frequency  peak  associated  with  the  smaller  squares,  as  expected.  Additionally,  the  higher 
frequency  peak  has  a  consistently  larger  magnitude,  most  likely  due  to  the  higher  density 
of  10  pm  squares.  It  is  also  apparent  that  the  strength  of  the  peak  at  6.5  TFIz  drops  faster 
compared  to  the  one  at  4  THz  as  the  thickness  of  the  dielectric  layer  is  increased.  The 
6.5  THz  peak  seems  to  be  much  more  sensitive  to  changes  in  dielectric  thickness 
compared  to  single-absorption  peak  configurations,  such  as  in  Figure  10(c). 


(a)  (b) 

Figure  14.  Unit  cell  (a)  and  absorption  characteristics  (b)  of  a  two  peak  metamaterial 
layer  utilizing  a  tile  configuration  to  maximize  square  density.  From  [23]. 


The  number  of  absorption  peaks  can  be  extended  from  dual-band  absorbers  by 

introducing  additional  square  sizes.  The  geometry  in  Fig.  15  (a)  has  three  different 

squares  (11,  14,  and  17  pm)  and  maintains  symmetry  in  two  directions,  thus  remaining 

polarization-independent  at  normal  incidence,  while  also  allowing  for  one  of  the  squares 

(in  this  case  the  14  pm  squares)  to  have  a  20%  higher  density,  enhancing  absorption  at 

one  of  the  corresponding  resonant  frequencies.  This  enhancement  can  be  useful  to 

increase  absorption  selectively  if  desired.  As  with  the  dual-band  configuration,  three 

distinct  peaks  associated  with  each  square  size  can  be  resolved  in  the  measured 

absorption  spectrum  in  Fig.  15(b).  In  this  case,  the  middle  peak  at  5  THz  has  a  higher 
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square  density  associated  with  it  and  an  absorption  peak  of  similar  magnitude  compared 
to  the  one  corresponding  to  the  larger  squares  at  4  THz.  The  overall  absorption  was  found 
to  be  lower  than  that  obtained  for  single-band  and  dual-band  metamaterials,  likely  due  to 
the  lower  density  of  particular  square  sizes  in  this  configuration. 


Frequency  (THz) 

(b) 


Figure  15.  Unit  cell  (a)  and  measured  absorption  (b)  for  a  triple  band  metamaterial 
configuration  with  square  sizes  of  1 1,  14,  and  17  pm  showing  three  clear 

peaks.  From  [23]. 


As  metamaterial  absorbers  have  a  narrowband  absorption  peak,  Kirchoff  s  Law 
demands  that  they  also  possess  a  narrowband  thermal  emission  peak  [50],  [52].  Indeed, 
the  emissivity  of  the  metamaterials  should  follow  their  absorptivity  exactly  for  thermal 
equilibrium  to  be  upheld.  As  described  earlier,  this  was  tested  using  the  metamaterial 
absorbers  as  an  external  source  for  the  FTIR  and  comparing  it  to  a  blackbody.  Samples 
A,  B,  and  C  consisted  of  single  band  absorbers  with  17  pm,  13  pm,  and  9  pm  squares, 
respectively.  All  samples  had  21  pm  periodicity  and  1.6  pm  thick  dielectric  spacer  of 
stoichiometric  SiC>2.  The  emission  spectrum  of  a  reference  emitter  (a  silicon  wafer  coated 
with  black  carbon)  was  used  to  extract  the  emissivity  of  the  metamaterial  samples  as  a 
function  of  THz  frequency.  The  emissivity  spectra,  shown  in  Figure  16  were  obtained  by 
dividing  the  sample  data  by  the  corresponding  data  obtained  using  the  blackened  emitter 
heated  to  the  same  temperature  (400  °C). 
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Figure  16.  Measured  emissivity  of  the  metamaterial  samples  A,  B  and  C  at  400  °C. 

Emissivity  exhibits  peaks  at  4.1,  5.4  and  7.8  THz,  respectively.  From  [24], 


It  is  clear  that  the  emission  spectra  in  Figure  16  agree  well  with  the  measured 
absorption  data  in  Figure  16  as  expected  from  Kirchhoff  s  law.  The  emission  bandwidth 
at  half  maximum  for  each  metamaterial  structure  is  approximately  1  THz. 
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The  total  power  emitted  into  each  THz  band  per  unit  area  (P)  can  be  estimated  by 
integrating  Planck’s  blackbody  radiation  formula  (P(v))  multiplied  by  the  measured 
emissivity  (s(v))  within  the  emission  bandwidth  as: 


P(v) 


A\ 

c  J 


:(v), 


exp 


hv/ 


kT  _ 


-dv 


-1 


(3) 


where  h  is  Planck’s  constant,  c  is  the  speed  of  light  in  vacuum  and  k  is  Boltzmann’s 
constant.  The  estimated  power  per  unit  area,  based  on  Eq.  (3),  for  samples  A,  B  and  C 
operated  at  400°C  is  11,  18  and  36  W/m2,  respectively. 


Figure  17.  Spectral  absorptance  of  the  metamaterial  samples  A,  B  and  C.  Solid  lines 
represent  the  measured  absorptance  with  peaks  at  4. 1,  5.4  and  7.8  THz, 
respectively.  Dashed  lines  represent  the  finite  element  simulation  results. 

From  [24]. 
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III.  METAMATERIALS  INTEGRATED  INTO  THERMAL 

SENSORS 


Note:  The  first  section  is  extended  from  a  paper  published  in  Optical  Materials 
Express  [25].  The  remaining  material  in  this  chapter  is  contained  entirely  within  a  paper 
published  in  Optics  Express  [26],  and  was  also  presented  at  the  SPIE  Photonics  West 
conference  in  February  2013. 

A.  THERMAL  SENSOR  PROPERTIES 

To  achieve  an  increase  in  performance  for  thermal  sensors,  we  must  first  consider 
their  performance  characteristics.  The  two  most  significant  parameters  are  the  thermal 
time  constant  (x)  and  responsivity  (Rth).  All  thermal  sensors  have  a  thermally  isolated 
sensitive  region  which  disperses  heat  to  the  environment  through  radiation,  convection 
and  conduction  as  illustrated  in  Figure  18. 


Sensor 

C 

G|eg 


Heat  sink 

P 

Figure  18.  Diagram  showing  heat  exchange  between  a  thermal  sensor  and  its 

environment.  Incoming  power  P  creates  a  gradient  of  AT  between  the  sensor 
and  the  heat  sink  at  temperature  T  in  addition  to  radiation  losses.  Convection 
losses  are  ignored  as  most  sensors  operate  at  low  pressures.  After  [20]. 


Radiation  losses  are  unavoidable  and  intrinsic  to  all  thermal  sensors  as  Kirchoff  s 

law  [50]  demands  that  any  object  that  absorbs  radiation  must  possess  the  same  emissivity, 

and  it  is  clear  that  THz  metamaterial  absorbers  are  no  exception  [24],  Convection  can  be 

controlled  through  the  gas  pressure  of  the  surrounding  environment  [52]  and  conduction 

through  materials  and  structural  design  of  the  sensors.  If  the  thermal  losses  can  be 

approximated  as  linear  with  respect  to  the  temperature  gradient  (AT)  between  the 
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thermally  isolated  sensitive  area  and  the  surrounding  environment,  they  can  be  lumped 
into  a  single  constant  tenn  called  the  thermal  conductance  ( Gt /,).  Therefore,  the  heat 
exchange  equation  for  a  thermal  sensor  absorbing  radiation  of  power  (Pt/,)  will  be: 


P*(t)  =  GAAT  +  Cm—-  (4) 

at 

where  Cth  is  the  thermal  capacitance  of  the  thermally  isolated  sensitive  area.  Solving  this 
differential  equation,  we  can  then  arrive  at  a  formula  for  the  responsivity  of  a  detector 
under  sinusoidal  illumination  power  [10]: 


R,h  = 


AT 


n 


Gj\+OfP 


(5) 


where  ;/  is  the  absorbance,  the  fraction  of  incident  power  (Po)  absorbed  by  the  sensitive 
region  of  the  detector,  r  is  the  thermal  time  constant  (r  =  Cth/Gth  )  and  co  is  the  angular 


modulation  frequency  of  the  incident  radiation.  The  graph  of  Rth  vs  Gth  in  Figure  19 
shows  that  reducing  thennal  conductance  is  always  advantageous  in  terms  of 
responsivity. 


Figure  19.  Responsivity  (Rth  [K/W])  vs  thermal  conductivity  (G,/,  [W/K])  in  a  thermal 
sensor  under  sinusoidal  illumination  power.  Clearly,  regardless  of  Gti„  R,h 
increases  with  higher  absorption  ( rj ),  lower  frequencies  {co),  and  smaller  C,i, 


32 


Indeed,  if  thennal  conductance  were  zero  the  sensor  would  simply  heat  up 
indefinitely  from  incident  radiation.  However,  real  sensors  do  not  have  an  infinite 
dynamic  range  and,  of  course,  Kirchoff  s  law  ensures  that  radiation  losses  cannot  be 
avoided  [50],  Most  thermal  sensors  are  designed  such  that  the  unavoidable  radiation 
losses  dominate  the  thermal  conductance  to  maximize  sensitivity  [10].  Therefore, 
integrating  metamaterials  will  primarily  influence  r  and  rj.  Clearly,  metamaterials  provide 
rj  close  to  100  %,  but  x  must  be  minimized  for  real  time  imaging  applications,  as  thermal 
detectors  typically  operate  in  a  frequency  range  where  cox  «  1  to  achieve  the  highest 
responsivity.  Therefore,  if  r  is  large,  to  must  be  small  making  real  time  imaging 
problematic.  This  means  that  because  Gth  is  already  minimized,  to  minimize  r  we  must 
minimize  Cth-  The  volume  of  the  materials  and  their  thermal  capacitance  will  determine 
the  total  Cth •  The  area  of  the  thermal  sensor  will  be  determined  by  the  application  and 
wavelengths  of  interest;  indeed,  more  sensor  area  can  result  in  a  higher  signal  from 
capturing  more  incident  power.  The  thickness  of  the  thermally  sensitive  area,  however, 
can  be  controlled  through  design  and  should  thus  be  as  small  as  possible.  This  is  the 
primary  reason  for  integrating  metamaterials  into  these  sensors,  the  structural  materials 
are  poor  THz  absorbers  by  themselves  [21],  [22],  and  excessively  thick  structures  will 
perform  poorly  in  a  real  time  imaging  setup. 

B.  BI-MATERIAL  SENSOR  CHARACTERISTICS 

As  the  previous  section  illustrates,  the  basic  physics  of  integrating  metamaterials 
into  thermal  sensors  is  similar  in  that  thermal  capacitance  and  absorption  are  the  primary 
concerns.  Of  course,  each  particular  sensor  will  have  its  own  challenges,  and  the  specifics 
of  these  challenges  for  bolometers  will  be  explored  in  the  following  chapter.  Initially,  we 
chose  bi-material  sensors  to  test  functionality  metamaterials  due  to  the  experience  of  my 
colleagues  in  their  fabrication  and  the  lab  equipment  available. 
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Like  all  thermal  detectors,  bi-material  sensors  respond  to  a  change  in  temperature 
of  their  absorbing  element;  in  this  case,  the  pixel  experiences  angular  deflection  per  unit 
power  ( dQ/dP )  and  this  responsivity  can  be  described  as  [26]: 

RJ±=R/A=—fi, — M  (6) 

dP  dT  gJ]  +  ott2  dT 

Clearly,  the  expression  is  simply  the  temperature  response  of  the  sensor  multiplied  by  the 
deflection  per  unit  temperature  change  (thermomechanical  sensitivity).  As  described 
before,  the  thermal  time  constant  limits  the  speed  of  the  sensor.  Noise  must  also  be 
examined  when  evaluating  performance;  in  bi-material  sensors  noise  arises  from  several 
different  sources  such  as  temperature  fluctuations,  background  fluctuations,  thermo¬ 
mechanical  resonances,  illumination  source  fluctuations  and  the  readout  system  [53].  The 
first  four  manifest  as  fluctuations  in  the  overall  sensor  deflection,  while  the  readout  noise 
depends  on  the  probing  mechanism.  In  a  practical  sense,  the  total  noise  of  the  complete 
detection  system  can  be  described  by  the  noise  equivalent  power  (NEP).  For  bi-material 
sensors,  NEP  can  be  defined  as  the  incident  radiant  power  that  produces  an  angular 
deflection  equal  to  detector’s  root  mean  square  (rms)  noise  [54],  [55], 

Fundamentally,  there  are  two  main  choices  when  designing  a  bi-material  sensor: 

materials  and  configuration.  Materials  should  be  fabrication-friendly,  exhibit  low  residual 

stress,  have  very  different  thermal  expansion  coefficients  and  strong  THz  absorption. 

Configurations  should  have  a  large  absorption  area,  good  thermal  isolation  to  increase 

sensitivity,  provide  a  reflective  surface  for  optical  readout,  and  have  a  small  thennal 

capacitance.  All  of  these  requirements  are  intrinsically  interdependent  making  the 

optimization  of  the  final  sensor  highly  dependent  on  the  intended  application. 

Nonetheless,  the  quest  to  achieve  high  perfonnance  THz  bi-material  detectors  starts  with 

dO/dT ,  defined  by  the  bimetallic  effect,  and  q,  which  is  maximized  by  using 

metamaterial  based  absorbers.  The  sensors  here  were  not  the  first  attempted,  in  our 

previous  work  [31]  a  complex  configuration  of  two  fourfold  legs  metalized  in  opposite 

sides  were  used.  Although  the  total  length,  and  thus  bi-material  effect,  was  higher  than 

that  of  the  current  design  (see  Figure  20(a)),  the  THz  sensitive  area  was  half  the  size  to 

keep  the  overall  pixel  sixe  around  200  x  200  pm2.  The  metallization  of  the  front  and  back 
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sides  of  alternate  legs  caused  the  structure  to  be  highly  stressed  after  release.  As  a  result, 
the  sensors  were  bent  as  much  as  30  degrees  out  of  the  substrate  plane,  which  is 
problematic  for  the  optical  readout.  The  configuration  reported  here  exhibits  several 
advantages,  which  are  highlighted  in  the  next  sections. 


Figure  20.  Bi-material  sensor,  (a)  3D  view  of  the  THz  bi-material  sensor  with  integrated 
metamaterial  absorber,  fabricated  on  a  Si  substrate,  (b)  Close  up  of  an 
isolated  bi-material  beam  which  length  is  h,  and  metal  and  dielectric 
thickness  fi  and  t2,  respectively.  Azieg  and  A9  are  the  linear  and  angular 
deflection  of  the  beam,  respectively.  From  [26]. 


To  increase  sensitivity,  it  is  important  to  optimize  the  bi-material  layer  thickness 
to  maximize  the  deflection  under  increasing  temperature.  The  linear  displacement  (Azieg) 
of  the  free  tip  of  a  bi-material  beam,  as  depicted  in  Fig.  19(b),  was  first  quantified  by 
Timoshenko  [56]  and  further  explored  by  other  groups  [53],  [57].  If  the  linear 
displacement  is  much  smaller  than  the  length  of  the  bi-material  beam  ( lb ),  which  is  true 

in  our  case,  the  angular  deflection  due  to  temperature  change  ( dQ/dT )  or 
thermomechanical  sensitivity  can  be  estimated  using: 
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where  t  represents  thickness,  a  is  the  thermal  expansion  coefficient  and  E  is  the  Young’s 
modulus.  The  indices  1  and  2  are  used  to  represent  materials  1  and  2,  respectively.  Table 
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1  lists  some  of  the  most  common  MEMS  materials  along  with  their  structural,  thermal 
and  electrical  characteristics. 


When  the  bi-material  legs  are  connected  to  a  freestanding  flat  absorber  (see 
Figure  21(a)),  the  sensor  angular  deflection  is  approximately  equal  to  AO  (see  Figure 
20(b)).  The  effect  can  be  further  amplified  by  adding  multifold  legs  with  alternate  bi¬ 
material  segments  [31].  However,  such  a  configuration  also  magnifies  the  bending  due  to 
residual  stress  after  release. 


Material 

Si 

SiNx 

Si02 

A1 

Au 

Young ’s  Modulus 
E(xl06  Pa) 
Expansion 

100 

180 

68 

70 

77 

Coefficient 
a  (x  10-  K1) 
Thermal 

2.7 

2.1 

0.4 

25 

14.2 

Conductivity 
g  (Wm'/C1) 

130 

19 

1.4 

237 

296 

Heat  Capacity 
c  (J  kg'1  K1) 

750 

691 

703 

900 

129 

Density 

p  (  x  10-  kg  m  3) 
Electric 

2330 

2400 

2200 

2700 

19300 

Conductivity 
a  (  x  106  S  m1) 
THz 

10 

37 

refractive  index 

* 

n 

3.48-0.  li 

2. 1— 0.025i 

2.0-0.2i 

Table  1.  Properties  of  standard  MEMS  materials [5 8],  [59].  From  [26]. 


It  is  clear  from  Figure  21  that  the  Al/SiCk  combination  produces  the  highest 
sensitivity  with  the  maximum  occurring  when  the  metal  thickness  is  approximately  one- 
half  of  the  dielectric  thickness.  Non-stoichiometric  SiNx  can  provide  less  stressed  layers 
than  SiCF  [60],  however,  silicon-rich  SiOx  can  be  deposited  with  much  lower  stress  than 
SiC>2,  while  preserving  most  of  the  thermomechanical  and  electro-optical  properties. 
During  sensor  fabrication,  testing  layers  of  non-stoichiometric  SiOx  and  stochiometric 
SiC>2  layers  with  the  same  thickness  were  deposited  on  Si  substrates  with  intrinsic  stress 
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on  the  order  of -13  MPa  and  -140  MPa,  respectively.  By  selecting  SiOx  and  A1  (both 
standard  MEMS  materials),  it  is  possible  to  maximize  dO/dT  while  simultaneously 
alleviating  some  of  the  excessive  residual  stress  related  deformation  observed  in  the 
sensors  fabricated  in  [31],  Furthermore,  SiOx  and  A1  exhibit  electro-optical  properties 
that  are  suitable  for  highly  efficient  metamaterial  absorbers,  as  shown  in  the  previous 
chapter. 


—  Al/Si02  -- 

—  Au/Si02 

- AI/SiN,  -- 

--- Au/SiN 

- Al/Si 

--- Au/Si 

•  AI/SiOx  (FE/ Experimental) 

o  dr 


_Bi-material  Leg 


Az„ 


A8 


Az, 


leg 


(a) 


(b) 


Figure  2 1 .  Thermomechanical  deflection  of  the  bi-material  sensor,  (a)  Freestanding  flat 
THz  absorber  connected  to  a  bi-material  beam,  whose  length  is  h,  and  metal 
and  dielectric  thickness  fi  and  t2, respectively.  Azabs  is  the  total  linear 
displacement  and  AQ  is  the  angular  deflection  of  the  absorber,  (b) 
Thennomechanical  sensitivity  ( d6/dT)  of  the  structure  of  part  (a),  calculated 
using  Equation  (4)  for  all  combinations  of  metal/dielectric  of  Table  1  where 
ti  varies  from  10  to  800  mn  and  t2  is  fixed  in  1.1  pm.  The  circular  marker 
shows  FE  and  experimental  results  for  ti  =  170  mn.  Solid  lines  are  for  A1 
combinations  while  dashed  lines  are  for  Au.  From  [26]. 


C.  BI-MATERIAL  THZ  SENSOR  DESIGN 

Bi-material  terahertz  detectors  were  designed  using  a  metamaterial  structure 
optimized  to  absorb  at  3.8  TFIz  radiation  form  a  QCL  available  to  us.  Relatively  large 
pixel  dimensions  were  chosen  to  increase  the  absorption  area  and  simplify  the  fabrication 
and  characterization  process.  Thennal  conductance  was  intentionally  varied  among  the 
designs  while  thermal  capacitance  remained  essentially  constant  (see  Table  2).  Figure  22 
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shows  the  geometrical  details  of  three  sensors  (A,  B,  and  C)  with  different  thermal 
conductances. 
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Figure  22.  Structural  parameters  of  the  three  bi-material  THz  sensors  in  pm.  (a)  Top  view 
of  sensor  A  showing  all  dimensions,  (b)  Top  view  of  sensors  B  and  C, 
showing  the  differences  in  sizes  of  the  thermal  insulator  anchors,  (c)  Vertical 
cut  of  the  sensor  structure.  From  [26]. 


The  detectors  consist  of  a  square  metamaterial  sensing  element  in  the  center 
connected  to  two  symmetrically  located  rectangular  bi-material  legs.  The  entire  structure 
is  then  connected  to  and  thermally  isolated  from  the  substrate  by  folded  SiOx  anchors 
with  varied  dimensions  as  shown  in  Figure  22.  The  thickness  of  the  A1  ground  plane  and 
squares  is  100  nm  while  the  bimaterial  legs  have  a  170  nm  layer  of  A1  on  the  top  side. 
The  structural  SiOx  is  1 . 1  pm  thick  everywhere.  The  thermal  conductance  (G)  of  all  the 
sensors  was  estimated  using  the  expression: 


Q  _  8  th 

l 


(8) 


where  gth  is  the  thermal  conductivity,  Ac  is  the  cross-sectional  area  and  /  is  the  length. 

Since  the  dimensions  of  the  thermal  isolation  sections  are  different,  the  total  thermal 

conductance  was  estimated  by  adding  the  thermal  resistance  of  each  section.  The 

metalized  parts  are  considered  thermal  shorts  due  to  their  high  thermal  conductivity 

compared  to  that  of  SiOx.  The  thennal  conductance  via  radiation  loss  of  heat  is  found  to 

be  an  order  of  magnitude  lower  than  that  via  the  legs  due  to  low  emissivity  of  A1  and  the 
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THz  metamaterial  that  cover  most  of  the  sensor  surfaces.  Heat  dissipation  due  to 
convection  is  negligible  as  the  sensors  typically  operate  under  low  pressure  (in  a  vacuum 
chamber).  The  thennal  capacitance  was  estimated  using  the  expression: 

c  =  cthPAst  (9) 

where,  clh  is  the  material  thermal  capacity,  p  is  the  material  density,  As  is  the  surface 

area  and  t  is  the  layer  thickness.  The  thermal  capacitance  of  the  sensor  is  the  sum  of 
thennal  capacitances  of  the  SiOx  and  A1  layers.  The  material  parameters  used  for  the 
calculations  are  given  in  Table  1.  The  time  constant  (x  =  C/G)  was  also  estimated  for  each 
sensor  configuration  and  listed  in  Table  2  in  addition  to  other  parameters. 


Sensor 

Sensor  A 

Sensor B 

Sensor  C 

Property 

Anal 

FE 

Exp. 

Anal 

FE 

Exp. 

Anal 

FE 

Exp. 

Absorptance 

•1 

- 

0.96 

0.95 

- 

0.96 

0.95 

- 

0.96 

0.95 

Thermal 

Conductance 

1.6 

1.7 

2.2 

2.1 

9.3 

8.5 

G  (x  10-  WK1) 
Thermal 
Capacitance 

11.1 

12 

10.7 

12.5 

9.8 

11.9 

C  (x  10-  JK1) 
Time  Constant 

0.68 

0.7 

0.8 

0.47 

0.6 

0.5 

0.1 

0.14 

0.3 

t(s) 

Thermomechanic 
al  Sensitivity 

0.15 

0.19 

0.2 

0.15 

0.2 

0.2 

0.15 

0.25 

0.2 

dd/dT  (deg  K1) 
Responsivity 
dO/dP  (x  106  deg 

0.95 

1.1 

1.2 

0.65 

0.9 

0.8 

0.15 

0.25 

0.2 

W1) 

Noise  Equivalent 
Power 

(due  to  incident 

0.00 

s 

8.6 

0.00 

f. 

13 

0.01 

A 

45 

power) 

NEP(x  10-  W) 

o 

Table  2.  THz  bi-material  sensor  analytical,  numerical,  and  experimental  parameters. 

From  [26], 
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The  deformation  of  the  sensor  with  increasing  temperature  was  analyzed  using  the 
COMSOL  heat  transfer  module,  which  allows  a  uniformly  distributed  heat  flux  boundary 
to  be  placed  at  the  absorber  to  emulate  the  incoming  THz  power.  The  anchor  attachments 
to  the  substrate  are  fixed  and  set  at  a  constant  temperature  to  represent  the  heat  sink.  All 
other  boundaries  are  thermally  insulated  from  the  surroundings  and  free  to  move.  The 
program  computes  the  heat  transfer  equation  at  each  mesh  point  allowing  the  retrieval  of 
several  parameters,  such  as  temperature  distribution,  thennal  deformation,  etc.  For  steady 
state  simulations  the  total  incoming  heat  flux  was  conveniently  set  as  1  pW,  therefore  the 
thennal  deformation  and  temperature  distribution  can  be  directly  read  “per  unit  pW.”  The 
angular  deformation  can  be  directly  obtained  by  the  displacement  of  the  free  edges  of 
absorber  and  hence  d9/dT  can  be  estimated  using  the  temperature  difference  between  the 
absorber  and  heat  sink.  Also,  the  responsivity  (d6/dP)  of  the  sensors  can  be  obtained 
using  the  maximum  deformation  (steady  state)  and  the  incident  heat  flux  (1  pW). 
Furthermore,  thennal  conductance  can  be  estimated  using  Equation  8. 

Time  domain  simulations  were  performed  to  obtain  the  transient  response  of  the 
sensor  to  a  pulsed  heat  flux  allowing  the  retrieval  of  the  time  constant  of  the  sensors. 
Using  the  obtained  time  constant  and  thermal  conductance,  we  estimated  the  thermal 
capacitance  of  the  sensors.  The  calculated  and  simulated  parameters,  using  the  material 
properties  of  Table  1,  are  listed  in  Table  2  and,  in  general,  show  good  agreement.  Notice 
that  the  thermal  capacitance  values  obtained  by  FE  simulations  show  a  small  discrepancy 
as  they  increase  with  decreasing  sensor  mass.  This  is  most  likely  due  to  the  time  constant 
estimation,  which  is  more  susceptible  to  errors  as  it  decreases.  Figure  24(a)  and  24(b) 
show  the  deformation  plots  obtained  by  FE  simulation,  where  the  z-axes  are  scaled  up  20 
times  for  visual  purposes.  The  surface  color  scale  indicates  the  temperature  distribution 
and  it  is  the  same  for  all  sensors.  It  can  be  seen  in  Figure  24  that  sensor  A  deflects  more 
compared  to  sensors  B  and  C  under  the  same  incident  power  (1  pW)  primarily  due  to 
lower  thermal  conductance.  Figure  24(d)  shows  the  time  domain  simulations  where  the 
sensors  are  submitted  to  a  step  excitation  (black  solid  line)  of  1  pW  for  duration  of  8 
seconds.  The  vertical  axes  show  temperature  on  the  left  side  and  angular  deflection  on  the 
right  side. 
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Figure  23.  FE  simulations  showing  the  deformation  plots  of  sensor  A  (a),  B  (b)  and  C  (c) 
under  a  constant  1  pW  heat  flux.  The  z-axes  are  scaled  up  20  times  for  visual 
purposes.  The  surface  colors  indicate  the  temperature  distribution  according 
to  the  color  bar  on  the  left,  (d)  Time  domain  simulation  of  all  three  sensors 
under  a  1  pW  step  excitation  (black  line).  Temperature  change  and  angular 
displacement  are  shown  on  the  left  and  on  the  right,  respectively.  From  [26]. 


Noise  sources  intrinsic  to  the  detectors  were  also  considered  and  an  analysis 
similar  to  that  in  [53]  was  performed  to  determine  the  NEP.  The  expressions  given  by 
Equations  10  and  11  were  adapted  from  [53]  to  reflect  angular  deflection  fluctuations. 
The  primary  noise  sources  in  thermal  detectors  are  temperature  fluctuation,  background 
fluctuation  and  thermomechanical  noises.  The  spontaneous  fluctuation  in  angular 
deflection  (deg)  of  the  absorbers  caused  by  temperature  fluctuations  is  given  by: 

(10) 

\  /  ^ 

where  T  is  the  sensor  temperature,  ke  is  the  Boltzmann  constant,  G  is  the  total  thermal 
conductance  and  B  is  the  bandwidth,  which  can  be  set  to  unity.  The  background 
fluctuation  noise  can  be  obtained  by  replacing  the  total  thermal  conductance  in  Equation 
10  by  thermal  conductance  via  radiation  loss  of  heat.  However,  this  is  much  smaller  than 
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the  thermal  conductance  via  the  legs  and  its  contribution  to  noise  can  be  neglected.  The 
angular  deflection  (deg)  due  to  thermomechanical  noise,  knowing  that  the  detector 
operating  frequency  is  much  slower  than  the  mechanical  resonances  (few  kHz),  is  given 
by: 


360  l4k3GB 
7ilh  ]]  Qka>0 


(11) 


where  Q  is  the  quality  factor,  k  is  the  stiffness  and  too  is  the  resonant  angular  frequency  of 
the  mechanical  structure.  Using  the  eigenfrequecy  solver  in  the  COMSOL  structural 
mechanics  module,  the  first  resonant  frequency  and  stiffness  of  all  the  sensors  were 
estimated  and  found  to  have  values  3.5,  4.0  and  6.0  kHz  and  0.02,  0.025  and  0.04  Nm  1 
for  sensors  A,  B  and  C,  respectively.  Typical  Q  values  for  similar  structures  lie  between 
100  and  1000  in  vacuum  [53].  The  noise  was  estimated  and  as  expected,  the  dominant 
source  is  the  temperature  fluctuation  in  the  detector.  The  total  noise  intrinsic  to  the 
sensors  was  estimated  to  be  5.0,  4.0  and  2.0  pdeg.  The  NEP  values  of  the  three  sensors 
were  calculated  by  dividing  the  fluctuations  due  to  the  noise  by  their  respective 
responsivities,  and  are  listed  in  Table  2. 


D.  FABRICATION  AND  CHARACTERIZATION 

The  sensors  were  fabricated  using  standard  micromachining  technology.  First,  a 
100  nm  thick  aluminum  (Al)  film  was  deposited  on  a  300  pm  thick  silicon  (Si)  substrate 
by  e-beam  evaporation.  Then,  the  Al  layer  was  patterned  and  wet  etched  to  form  the 
absorber  ground  plane.  Next,  a  1.1  pm  thick  SiOx  layer  was  deposited  using  plasma 
enhanced  chemical  vapor  deposition  (PECVD)  at  300  °C,  followed  by  another  100  nm 
thick  Al  film.  The  second  Al  layer  was  then  patterned  and  plasma  etched  to  define  the 
absorber  metamaterial  squares.  Then  a  170  nm  thick  Al  layer  was  deposited,  patterned 
and  lifted  off  to  form  the  bi-material  legs.  The  sensor  structure  was  then  created  by 
reactive  ion  etching  of  the  SiOx  layer.  Finally,  the  structures  were  released  through 
backside  trenching  using  the  Bosch  etch  process.  Circular  openings  were  chosen  to 
ensure  release  of  the  structure  and  to  help  refine  the  Bosch  etch  recipe.  Figure  24(a)  and 
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24(c)  show  the  optical  profile  of  the  fabricated  sensor  A  along  with  micrographs  of  the 
three  sensors  and  an  array  of  sensors  of  type  A. 


Figure  24.  Fabricated  THz  bi-material  sensors,  (a)  3D  optical  profile  of  sensor  A  (the 
aspect  ratio  is  preserved),  (b)  Micrograph  of  an  array  of  sensor  A.  (c)  2D 
profile  taken  along  the  bimaterial  legs  direction  (y-profile)  with  the 
processing  direction  (z-profile)  scale  exaggerated  to  show  the  residual 
deflection  of  the  legs  (red  line)  and  absorber  (blue  line),  (c)  Micrographs 
showing  the  top  view  of  sensors  A,  B  and  C.  From  [26]. 


The  3D  profile  in  Figure  24(a)  is  an  actual  view  of  sensor  A  while  the  2D  profile 
in  Figure  24(c),  taken  along  the  bi-material  legs  direction  (y-profile),  is  exaggerated  in 
the  processing  direction  (z-profile)  to  show  the  residual  deflection  of  the  legs  and 
absorber  of  the  sensor.  The  measured  residual  deflection  of  the  absorber  is  approximately 
6°  for  the  sensors  A  and  B  and  8°  for  sensor  C.  It  is  easy  to  observe  that  the  absorber  is 
almost  flat  due  compensation  of  stresses  from  the  aluminum  layers  in  both  sides  of  the 
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SiOx  layer.  Due  to  the  deflection  of  the  sensors,  micrographs  shown  in  Figure  24(d)  are 
not  completely  focused  across  the  surface.  In  addition  to  the  sensors,  the  fabricated  wafer 
contains  an  area  of  10  x  10  mm  filled  with  the  same  metamaterial  structure  used  in  the 
sensors.  This  is  to  allow  accurate  measurement  of  the  absorption  characteristics  of 
metamaterial  used  in  the  sensors.  The  absorptance  of  the  metamaterial  film  was  measured 
as  described  in  chapter  II  and  compared  with  the  QCL  emission  characteristics  as  shown 
in  Figure  25(a).  A  good  match  between  the  absorptance  peak  position  of  the  metamaterial 
and  the  3.8  TFIz  QCL  emission  frequency  was  achieved  with  nearly  90%  absorptance. 
This  assured  that  the  sensors  absorbed  our  QCL  emission  with  high  efficiency. 


Figure  25.  (a)  Measurement  of  the  absorptance  spectra  of  the  THz  sensors  metamaterial 
structure  (blue  line)  compared  with  the  QCL  normalized  emission  (read  line), 
(b)  Measured  angular  deflection  (markers)  upon  temperature  chance.  Notice 
that  the  effect  is  linear  and  almost  indistinguishable  among  the  sensors, 
resulting  in  approximately  0.2  deg/K  thermomechanical  sensitivity.  From 

[26]. 


Next,  the  thermal  response  of  the  sensor  ( dd/dT)  was  measured.  The  temperature 
gradient  in  the  bi-material  section  of  the  leg  was  estimated  to  be  less  than  5%  of  that 
between  the  central  absorbing  element  and  the  substrate.  Thus,  the  bi-material  section  of 
the  leg  can  be  treated  as  thermally  shorted  allowing  the  measurement  of  the  thermal 
response  by  unifonnly  heating  the  sensor.  The  measurement  was  performed  by  attaching 
the  sensor  to  a  flat  resistive  heating  element  and  sweeping  the  temperature  from  303  to 
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313  K.  The  reflection  of  a  laser  diode  beam  from  the  backside  of  the  sensor’s  ground 
plane  was  projected  on  a  screen  and  the  angular  deflection  of  the  sensor  was  detennined. 
Angular  deflections  from  the  three  sensors  are  shown  in  Figure  25(b)  with  different 
markers.  The  deflections  are  almost  indistinguishable  because  the  detectors  have  the 
same  bi-material  leg  dimensions.  The  solid  line  is  a  linear  fit,  showing  that  the  response 
in  this  temperature  range  is  linear  and  approximately  0.2  deg/K,  which  is  slightly  higher 
than  the  estimated  values  (see  Table  2  and  Figure  21).  Figure  21  shows  that  the  thennal 
response  of  the  sensors  can  be  further  increased  by,  for  example,  decreasing  the  SiOx 
thickness  or  increasing  the  A1  thickness  on  the  legs.  Test  structures  fabricated  in  parallel 
with  these  sensors  showed  that  increasing  the  A1  thickness  on  the  legs  also  increases  the 
residual  stress.  Decreasing  the  dielectric  thickness  has  a  similar  effect  in  addition  to 
reducing  absorber  efficiency  [30].  Further  tuning  of  the  fabrication  process  is  necessary 
to  reduce  the  residual  stress  on  the  bi-material  legs  to  decrease  the  initial  bending  as 
depicted  in  Figure  7(a). 

Subsequently,  the  sensors  were  placed  in  a  vacuum  chamber  and  operated  at  a 
pressure  of  approximately  0.03  mTorr  to  minimize  the  heat  loss  by  convection  [52],  The 
QCL  was  kept  inside  a  cryostat  and  operated  at  around  15  K.  The  divergent  THz  beam 
passed  through  the  cryostat  Tsurupica  window  and  the  radiation  was  focused  by  a  40  mm 
polyethylene  lens  onto  the  sensors.  Both  Tsurupica  and  polyethylene  exhibit  reasonable 
transmission  (~65%)  in  the  THz  range.  The  QCL  was  operated  in  pulsed  mode  with  the 
pulse  width  fixed  at  5  ps  and  a  variable  pulse  rate  to  control  the  output  power.  The 
deflection  of  the  sensor  was  measured  using  the  same  procedure  described  earlier  for  a 
set  of  QCL  pulse  rates  ranging  from  nearly  zero  to  5  kHz.  The  absolute  power  that 
reaches  the  sensors  (incident  power)  is  estimated  using  the  responsivity  (dO/dP)  in 
Equation  6  along  with  the  calculated  thermal  conductance  and  measured  absorptance. 
Note  that  the  QCL  switching  frequency  and  duty  factor  must  be  taken  into  account  since 
the  sensors  can  only  respond  to  the  average  power.  Figure  26(a)  shows  the  measured 
angular  deflection  versus  incident  power  for  all  three  sensors  (colored  markers).  For  all  of 
the  sensors,  the  responsivity  values  estimated  analytically  tend  to  be  lower  than  that  of 
the  FE  and  experimental  values  mainly  due  to  the  underestimation  of  bimetallic  effect 
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(dd/dT),  by  Equation  7  (see  Figure  21).  As  expected  from  Equation  6,  responsivity  of  the 
sensors  was  found  to  decrease  with  increasing  thermal  conductance. 


Figure  26.  Responsivity  and  NEP  measurements,  (a)  Measured  angular  deflection  per 
varying  incident  power  for  all  three  sensors  (colored  markers)  Notice  that 
responsivity  increases  as  thermal  conductance  decreases,  (b)  Measured  output 
voltage  of  the  PSD  for  sensor  A  by  gating  the  QCL  output  at  200  mHz.  The 
power  incident  in  the  detector  is  shown  on  the  right  vertical  axis.  From  [26]. 


To  determine  NEP,  a  position-sensing  detector  (PSD)  was  added  to  the 
experimental  setup  to  read  the  deflection  at  low  power  levels.  The  NEP  was  then 
measured  for  each  detector  and  listed  in  Table  2.  Figure  26(b)  shows  measured  output 
voltage  of  the  PSD  for  sensor  A  by  gating  the  QCL  output  at  200  mHz.  It  is  important  to 
highlight  that  the  measurements  include  the  effects  of  QCL  power  fluctuations  and 
optical  readout  noise,  not  considered  in  the  theoretical  estimations  discussed  earlier.  The 
difference  between  the  measured  values  (3  orders  of  magnitude  higher)  and  the  estimated 
ones  (Table  2)  can  be  attributed  primarily  to  the  readout  noise.  The  QCL  power 
fluctuations  [61]  do  not  seem  to  contribute  to  the  observed  noise  since  the  noise  floor 
when  the  QCL  is  off,  shown  in  Figure  26(b),  is  similar  to  the  noise  observed  when  the 
QCL  is  on.  As  expected,  NEP  increases  from  sensor  A  to  C  due  to  decrease  in 
responsivity.  The  measured  NEP  values,  including  the  readout  noise  and  the  intrinsic 
noise  of  the  sensor,  can  be  translated  into  minimum  detectable  temperature  difference  on 
the  sensor,  found  to  be  approximately  50  mK  for  all  three  sensors.  This  value  is  similar  to 
those  of  bi-material  sensors  operating  in  the  IR  range  [53],  [58]. 
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The  time  domain  response  was  also  measured  using  the  PSD  and  the  results  for 
the  three  sensors  are  shown  in  Figure  27(a)  under  the  same  incident  power  with  the  QCL 
gated  at  500  mHz.  As  observed  in  Figure  26,  sensor  A  is  more  sensitive,  which  agrees 
with  the  predictions  and  previous  measurements.  Since  the  sensors  have  the  same  r\,  the 
same  absorbing  area,  same  materials,  the  same  dd/dT,  and  nearly  the  same  thermal 
capacitance,  speed  and  responsivity  are  completely  controlled  by  the  thermal 
conductance,  which  depends  on  the  anchor  geometry. 


Figure  27.  Time  and  frequency  domain  measurements,  (a)  Time  responses  of  sensors  A, 
B  and  C  measured  under  the  same  incident  power  with  the  QCL  gated  at  1 
Hz.  Noticed  that  sensor  A  is  more  sensitive  and  slower,  which  agrees  with  the 
predictions  and  previous  measurements,  (b)  Normalized  frequency  responses 
for  the  three  sensors  (colored  lines).  The  time  constants  were  retrieved  by 
taking  the  inverse  of  the  3  dB  frequencies  that  are  1.2,  2.1  and  3.2  rad/s  for 
sensor  A,  B  and  C,  respectively.  From  [26]. 


The  time  constant  of  the  sensors  was  determined  by  sweeping  the  QCL  gating 
frequency  from  50  mHz  to  30  Hz  and  recording  the  PSD  peak  to  peak  voltage.  The 
normalized  frequency  responses  for  the  three  sensors  are  shown  in  Figure  10(b).  The  time 
constants  were  retrieved  from  Figure  27(b)  by  taking  the  inverse  of  the  3  dB  frequency 
and  included  in  Table  2.  In  general,  the  measured  time  constants  agree  well  with  the  FE 
estimations,  while  the  analytical  approach  underestimates  this  parameter.  The  sensors  are, 
apparently,  slow  for  real  time  imaging;  however,  pixel  size  can  be  reduced  to  increase  the 

speed  of  operation.  A  tradeoff  between  speed  and  sensitivity  has  to  be  made  according  to 
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specific  applications.  Although  the  fabricated  sensor  arrays  do  not  have  high  spatial 
resolution,  their  imaging  capabilities  were  probed  by  a  CCD  camera  with  coaxial 
illumination  as  schematically  illustrated  in  Figure  28(a)  [62],  [63].  An  image  of  the  QCL 
beam,  gated  at  500  mHz,  is  shown  in  Figure  28(c)  compared  with  the  same  image 
obtained  using  a  commercial  IR  microbolometer  camera  with  THz  optics  [64]  (Figure 
27(b)).  The  focal  plane  array  of  the  IR  camera  has  30  pm  pitch  of  and  can  resolve  the 
rings  associated  with  the  QCL  beam.  Our  sensor  array,  on  the  other  hand,  has  a  430  pm 
pitch  and  cannot  resolve  the  rings;  nevertheless,  it  gives  a  raw  image  that  clearly  shows 
where  the  energy  is  concentrated  and  the  circular  shape  of  the  THz  beam.  Finally,  it  is 
possible  to  optimize  THz  bi-material  sensor  configuration,  size,  fabrication  processes  and 
readout  to  achieve  real  time  imaging. 
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(b)  (c) 


Figure  28.  QCL  beam  imaging,  (a)  Optical  readout  used  to  record  videos  [46]  and  the 
snap  shot  showed  in  part  (c).  The  images  were  recorded  using  background 
subtraction  to  suppress  to  the  effects  of  the  residual  stress  of  the  sensors,  (b) 
Image  obtained  using  a  30  pm  pitch  commercial  IR  microbolometer  camera 
with  THz  optics,  (c)  Image  (Media  1)  of  the  same  QCL  beam,  gated  at  500 
mHz,  obtained  using  an  array  of  sensor  A  with  430  pm  pitch  using  the 
readout  depicted  in  part  (a).  Notice  that  since  the  pitch  of  our  sensor  are  one 
order  of  magnitude  higher  than  the  IR  camera,  it  cannot  resolve  the  rings 
associated  with  the  QCL  beam,  showed  in  part  (b).  From  [26]. 
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IV.  INTEGRATING  METAMATERIALS  INTO 
MICROBOLOMETER  PIXELS 


Note:  Other  than  the  last  section,  this  chapter  is  published  practically  unaltered 
in  Optical  Materials  Express  [25]. 

A.  CHALLENGES  TO  BOLOMETER  INTEGRATION 

As  the  previous  chapters  have  described,  metamaterials  can  make  effective  thin 
THz  absorbers  for  thermal  sensors.  Bi-material  sensors  demonstrated  in  Chapter  III  the 
metamaterial  THz  absorber  and  bi-material  legs  are  physically  separated  and  do  not 
interfere  with  each  other’s  function.  In  the  case  of  microbolometers,  a  temperature 
sensitive  layer  must  be  included  in  the  central  absorbing  area,  or  microbridge  which  can 
affect  the  efficiency  of  the  metamaterial  absorber.  This  resistive  bolometric  layer  can  be 
made  out  of  several  different  materials,  such  as  titanium  (Ti)  and  vanadium  oxides 
(VaOx).  In  Chapter  III,  it  was  shown  that  the  absorbing  layers  should  be  as  thin  as 
possible  to  reduce  thennal  capacitance  (Cth),  and  thus  the  thinnest  choice  would  be  one  in 
which  the  bolometric  layer  is  simply  integrated  into  the  middle  of  the  dielectric  spacer  as 
in  Figure  29  or  serve  as  the  ground  plane  of  the  metamaterial.  Unfortunately,  the  latter 
option  leads  to  a  significant  reduction  in  absorption  as  THz  will  easily  penetrate  high 
resistivity  bolometric  materials  and  thus  defeats  the  purpose.  Therefore,  the  fonner  option 
was  investigated,  although  some  other  possible  configurations  will  be  examined. 

B.  FABRICATION  AND  MEASUREMENT 

The  approach  used  in  this  work  involves  a  metamaterial  structure  with  a 
bolometric  layer  embedded  in  the  dielectric  layer  as  schematically  illustrated  in  Figure 
29.  The  aim  of  the  work  in  this  chapter  is  to  detennine  the  optimal  location  of  the 
bolometric  layer  within  the  dielectric  to  achieve  the  highest  possible  THz  absorption. 
Three  wafers  containing  metamaterial  absorbers  with  Ti  as  the  bolometric  layer  located  in 
the  dielectric  spacer  at  varying  depths  along  with  one  without  a  bolometric  layer  were 
fabricated  using  standard  MEMS  microfabrication  techniques,  similar  to  those  used  in 
Chapter  II.  The  fabrication  sequence  includes  deposition  of  an  85  nm  thick  layer  of  A1  on 
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4”  Si  substrates  using  e-beam  evaporation.  Following  this,  a  layer  of  low  stress 
nonstochiometric  silicon-rich  silicon  oxide  (SiOx)  with  different  thicknesses  was 
deposited  using  plasma  enhanced  chemical  vapor  deposition  (PECVD). 


Figure  29.  Schematic  (a)  and  top  view  micrograph  (b)  of  a  metamaterial  absorber  with  a 
bolometric  layer  embedded  in  the  dielectric.  The  square  size  and  pitch  were 
selected  to  have  an  absorption  peak  within  the  1-10  THz  band.  From  [25]. 

Wafers  A-D  possessed  1270  mn,  765  nm,  297  mn,  and  1590  mn  of  SiOx, 
respectively,  confirmed  through  optical  interferometry.  Then,  a  Ti  layer  was  sputter 
deposited  using  Ar  plasma  onto  wafers  A-C  and  a  Si  test  wafer  with  a  SiNx  isolation  layer 
to  measure  conductivity.  The  measured  thickness  of  Ti  was  found  to  be  about  10  nm  and 
confirmed  to  have  a  sheet  resistance  of  about  220  D/sq,  which  corresponds  to 
conductivity  of  4.45  X  105  S/m.  Following  the  deposition  of  Ti,  an  additional  SiOx  layer 
with  thicknesses  of  297  nm,  765  nm,  and  1270  nm  were  deposited  using  PECVD  on 
wafers  A-C,  respectively.  This  resulted  in  nearly  equal  total  spacer  thicknesses,  including 
Ti  layer,  of  1577  nm,  1540  nm,  1577  nm,  and  1590  nm  for  wafers  A-D.  A  second  85  nm 
thick  layer  of  A1  was  then  deposited  onto  wafers  A-D  using  e-beam  evaporation  and  was 
then  patterned  using  photolithography  and  sputter-etched  with  Ar  plasma  to  form  an 
arrays  of  squares  as  in  Chapter  II,  with  the  absorption  characteristics  of  the  15  pm  side 
and  2 1  pm  pitch  configuration  shown  here. 
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As  described  in  Chapter  II,  square  size  mainly  influences  the  resonant  absorption 
frequency,  while  dielectric  thickness  and  periodicity  strongly  affects  absorption 
magnitude.  This  particular  configuration  was  chosen  for  its  resonant  frequency  in  the  1- 
10  THz  range,  and  other  configurations  with  varying  resonant  frequencies  can  be  found 
in  Chapter  II.  All  wafers  were  intentionally  over-etched  approximately  150  mn  into  the 
SiOx  to  assure  complete  insulation  between  square  elements.  The  composition  of  layers 
of  wafers  A-D  are  summarized  in  Table  3.  To  detennine  the  optical  properties  of  these 
metamaterial  absorbing  layers  in  the  THz  spectral  range,  reflection  measurements  were 
taken  using  a  Thermo-Nicolet  Nexus  870  Fourier  Transform  Infrared  Spectrometer 
(FTIR)  with  a  globar  source  and  PIKE  Technologies  MappIR  accessory.  An  Au-coated  Si 
wafer  was  used  to  establish  the  background  for  the  reflectance  measurements.  The  two  85 
mn  thick  A1  layers  block  any  significant  transmission  due  to  the  40  mn  skin  depth  at 
around  4.7  THz,  resulting  in  over  two  orders  of  magnitude  reduction  of  power  due  to 
exponential  attenuation.  The  periodicities  used  are  less  than  one  third  of  the  free  space 
wavelength,  resulting  in  negligible  higher  order  scattering  [43].  As  before,  this  allows 
A+T  +  R  =1  to  be  simplified  to  A  =1-  R  .  The  resulting  absorption  spectra  for  samples 
A-C  are  plotted  in  Figure  30(a). 


Wafer 

Top  SiOx 

Bottom  SiOx 

Total 
Spacer ‘ 

Max. 

Absorption1' 

A 

297 

1270 

1577 

62% 

B 

765 

765 

1540 

88% 

C 

1270 

297 

1577 

94% 

D 

— 

— 

1590 

98% 

"Includes  10  mn  of  Ti. 

^For  15  pm  squares  with  21  pm  periodicity 

Table  3.  Thickness  of  the  dielectric  spacer  layers  (in  mn)  of  each  sample.  From  [25]. 
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(a) 


(b) 


Figure  30.  (a)  Experimental  (solid)  and  FE  model  (dashed)  absorption  spectra  for 

samples  A  (black),  B  (red),  and  C  (green),  (b)  Theoretical  absorption  curves 
for  sample  A  for  a  set  of  sheet  resistivites  in  the  bolometric  layer.  Square  size 
is  15  pm  with  21  pm  periodicity  for  both  plots.  From  [25]. 


C.  RESULTS  AND  DISCUSSION 

It  is  clear  that  in  Figure  30(a),  as  the  Ti  conductive  layer  gets  closer  to  the 
patterned  top  A1  layer,  the  THz  absorption  of  the  metamaterials  is  significantly  reduced 
and  slightly  red  shifted.  This  is  expected  conceptually  as  the  Ti  layer  would  function  as  a 
poor  ground  plane  if  the  lower  SiOx  and  A1  layers  were  removed.  The  redshift  in  the 
resonant  absorption  frequency  is  larger  than  the  frequency  shifts  observed  in  Chapter  II 
for  thinner  dielectric  thicknesses,  suggesting  that  a  simple  reduction  of  dielectric 
thicknesses  will  not  explain  these  absorption  spectra.  To  model  the  THz  absorption 
characteristics  of  the  metamaterial  film  with  an  embedded  bolometric  layer,  finite 
element  (FE)  simulations  were  performed  using  the  RF  module  of  COMSOL 
multiphysics  software  explained  in  greater  detail  in  Chapter  II.  The  periodic  geometry  of 
the  metamaterial  absorbing  film  allows  for  a  single  unit  cell  to  be  modeled  with  periodic 
boundary  conditions.  Internal  ports  or  scattering  boundary  conditions  allow  incoming 
plane  waves  at  specific  frequencies  to  enter  the  model.  The  THz  absorption  can  be 
calculated  by  subtracting  off  reflection  and  transmission  coefficients  calculated  from  the 
S  parameters.  Alternatively,  resistive  heating  can  be  used  to  calculate  absorption.  This 
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method  is  especially  useful  for  calculating  the  absorption  of  elements  that  are  not 
thermally  linked,  such  as  a  ground  plane  disconnected  from  a  bolometer  microbridge. 
The  aluminum  layers  were  modeled  using  a  conductivity  of  1  x  107  S/m,  confirmed 
through  four-point  probe  and  stylus  profilometer  measurements.  SiOx  was  modeled  as  a 
dielectric  with  a  complex  index  of  refraction  of  1.95  +  0.025i,  adjusted  from  the  SiC>2 
values  used  in  Chapter  II  based  on  the  experimental  results  therein.  Ti  was  modeled  as 
having  a  sheet  resistance  of  220  £2/sq  based  on  four-point  probe  measurements.  The  exact 
thickness  of  the  Ti  is  irrelevant  as  long  as  the  sheet  resistance  is  well  known  [48]  and  was 
estimated  at  10  mn  from  previous  depositions.  As  in  Chapter  II,  the  over-etching  of  SiOx 
is  included  in  the  models,  producing  a  slight  blue  shift  in  the  resonant  absorption 
frequency  due  to  an  effectively  smaller  index  of  refraction.  Figure  30(a)  shows  the 
comparison  of  simulated  and  measured  data  for  the  three  samples.  The  simulated  data  is 
able  to  capture  the  effect  of  the  location  of  Ti  layer  on  the  absorption  reasonably  well.  E- 
field  magnitude  plots  for  these  simulations  are  shown  in  Figure  31,  all  using  the  same 
color  scale  for  E-field  intensity. 


(b) 


(c) 


Figure  3 1 .  E-field  magnitude  plots  from  COMSOL  FE  simulations  of  1 5  pm  squares 
with  21  pm  pitch  at  4.7  THz  for  samples  A  (a),  B  (b),  and  C  (c).  All  plots  use 
the  same  color  scale,  with  warmer  colors  indicating  more  intense  E-fields. 
Material  types  of  different  regions  are  labeled  in  (a).  E-field,  H-field,  and 
propagation  vectors  of  the  incident  THz  waves  are  shown  on  (b).  From  [25]. 
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It  is  clear  that  the  E-fields  do  penetrate  the  Ti  layer  but  not  the  ground  plane  in  all 
cases,  and  that  the  Ti  is  indeed  acting  as  a  poor  ground  plane.  The  E-fields  also  become 
less  intense  as  the  Ti  layer  nears  the  square  array,  indicating  that  it  is  negatively  affecting 
the  resonance.  In  addition,  absorption  of  sample  A  is  also  simulated  for  a  set  of 
conductivities  of  the  bolometric  layer  as  shown  in  Figure  30(b).  It  can  be  seen  that,  if  the 
conductivity  of  this  layer  is  reduced,  the  influence  of  the  layer  on  the  absorption  of  the 
metamaterial  decreases,  as  expected.  Since  typical  high  TCR  bolometric  layers  have 
sheet  resistances  in  the  kfl  range,  we  should  expect  that  a  spacer  of  300  nm  of  SiOx 
between  the  bolometric  layer  and  the  patterned  A1  layer  will  be  sufficient  for  achieving 
high  THz  absorption.  The  results  indicate  that  the  two  most  important  parameters  that 
affect  the  absorption  peak  height  are  the  conductivity  and  position  of  the  bolometric  layer 
between  the  patterned  A1  layer  and  ground  plane  of  the  metamaterial.  The  data  in  Figure 
30  also  seems  to  show  that  within  the  range  of  parameters  used,  there  is  some 
equivalence  in  the  effects  of  sheet  resistance  and  position  of  a  thin  bolometric  layer  on 
the  absorption  spectrum.  Qualitatively,  less  conductive  layers  should  allow  for  greater 
penetration  of  E-fields  and  thus  more  of  the  desirable  resonance  between  the  square  array 
and  the  ground  plane.  Also,  a  high  conductivity  layer  near  the  ground  plane  or  a  low 
conductivity  layer  near  the  square  array  should  both  converge  to  the  original  absorption 
spectrum. 

D.  OTHER  POSSIBLE  CONFIGURATIONS 

Now  that  we  have  established  that  the  impact  on  absorption  is  reduced  as  the 
sheet  conductivity  is  reduced  and  as  it  is  moved  closer  to  the  ground  plane,  we  can 
evaluate  other  possible  configurations  investigated  by  other  researchers  as  depicted  in 
Figure  32.  Maier  et  al.  fabricated  a  structure  analogous  to  Figure  32(a)  [65]  and  Figure 
32(c)  [32]  for  the  mid-IR  range  (3-8  pm)  range.  They  note  that  the  bolometric  layer,  in 
their  case  polysilicon  (poly-Si),  interfered  with  their  metamaterial  significantly  for  the 
configuration  in  Figure  32(c).  Certainly,  the  configuration  in  Figure  32(a)  eliminates 
interference  between  metamaterial  and  the  bolometric  layer  since  they  are  optically 
separated  since  the  ground  plane  pennitting  no  significant  THz  transmission.  However, 

when  comparing  THz  metamaterials  to  mid-IR  metamaterials,  an  important  distinction 
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emerges.  Note  that  in  Chapter  II,  the  optimal  dielectric  thickness  for  high  THz  absorption 
was  approximately  1  pm,  while  in  the  case  of  Maier  et  al.  the  optimal  dielectric  thickness 
can  be  as  thin  as  250  nm  [65]  also  seen  in  other  similar  work  [42].  Thus,  the 
configuration  in  Figure  32(a)  is  not  suitable  for  THz  applications  due  to  increase  of 
thermal  capacitance  when  extra  layers  added  for  separating  the  bolometric  layer. 


(a)  (b)  (c) 


Figure  32.  Three  possible  configurations  for  integrating  a  conductive  layer  into  a  THz 
absorbing  metamaterial  using  a  metal,  bolometric  layer,  and  a  dielectric  with 
the  bolometric  layer  below  (a)  and  within  (c)  the  metamaterial  layer. 
Configuration  (c)  thermally  disconnects  the  ground  plane.  The  configuration 
we  used  is  shown  in  (b)  with  Ti. 

Of  course,  thickness  can  be  reduced  with  some  loss  in  absorption  if  required  and 
controlled  over-etching  can  eliminate  some  of  the  dielectric  spacer.  However,  these 
methods  will  still  reduce  the  performance  of  the  microbolometer  pixel;  this  is  why  the 
configuration  in  Figure  32(b)  was  emphasized  for  THz  metamaterials  in  addition  to  the 
fact  that,  as  with  the  bi-material  pixels,  the  metal  layers  on  both  sides  of  the  dielectric 
will  compensate  for  each  other  in  terms  of  residual  stress.  The  structure  depicted  in 
Figure  32(c)  does  present  an  interesting  opportunity,  albeit  with  some  fabrication 
challenges.  In  this  configuration,  the  thickness  of  the  metamaterial  absorbing  layer 
includes  the  air  gap  thus  allowing  the  thermal  capacitance  of  the  microbridge  to  be 
relatively  small.  The  issue  here  is  shown  in  Figure  5  in  Chapter  II,  the  FE  model  indicates 
that  the  ground  plane  has  a  non-negligible  power  dissipation  of  over  20%  of  the  absorbed 
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THz  at  the  resonant  frequency.  As  the  dissipation  in  any  material  will  depend  on  the 
electric  field  and  conductivity  or  complex  refractive  index  of  that  material  (assuming 
magnetic  losses  are  negligible),  it  should  be  possible  to  effectively  control  where  the 
absorption  occurs.  Indeed,  if  the  ground  plane  were  constructed  out  of  a  material  with  a 
high  THz  conductivity,  dissipation  in  that  layer  would  be  less  than  the  other  layers, 
minimizing  this  problem.  Therefore,  while  the  structure  in  Figure  32(b)  is  likely  the 
safest  configuration  to  fabricate  microbolometers  with  high  perfonnance,  Figure  32(c) 
would  allow  for  an  even  smaller  thennal  capacitance. 
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V.  CONCLUSION 


The  research  presented  here  demonstrates  the  ability  of  metamaterials  to  enhance 
thennal  sensors  operating  in  the  THz  region  of  the  spectrum  by  improving  their 
absorption  characteristics.  Finite  element  models  provided  key  insight  leading  to 
fabricated  metamaterial  absorbers  and  bi-material  sensors  that  performed  well  in  the  3-8 
THz  range. 

The  THz  metamaterial  absorbers,  composed  of  a  metal  layer  patterned  into  a 
periodic  array  of  squares  separated  from  a  ground  plane  of  continuous  metal  by  a 
dielectric  spacer  were  found  to  have  nearly  100%  absorption  at  a  specific  THz  frequency. 
The  peak  absorption  frequency  was  inversely  proportional  to  the  dimension  of  the 
squares.  In  addition,  metamaterial  THz  absorbing  layers  with  one,  two,  and  three 
absorption  peaks  we  fabricated  using  complex  unit  cell  structures  instead  of  single  square 
unit  cells.  The  use  of  standard  microfabrication  materials,  A1  and  SiOx  for  the  fabrication 
of  metamaterials,  allows  them  to  be  incorporated  into  the  construction  of  thermal  sensors 
such  as  microbolometers  or  bi-material  sensors.  The  results  indicate  that  these 
metamaterial  absorbers  can  be  used  as  absorbing  layers  in  thermal  imagers  paired  with 
narrow-band  sources  to  significantly  enhance  their  sensitivity  through  resonant 
absorption  at  the  emitter  frequency.  In  addition,  dual-  and  triple-band  configurations  can 
be  used  with  multiple  lasers  or  broader  sources,  as  desired. 

The  metamaterials  were  successfully  integrated  into  bi-material  pixels.  Sensor 
materials  and  configurations  were  chosen  in  order  to  maximize  responsivity.  The 
combination  of  favorable  thermal,  mechanical  and  optical  properties  of  the  MEMS 
fabrication-friendly  materials  SiOx  and  A1  were  advantageous.  Analytical  and  FE  models 
were  used  to  predict  the  perfonnance  of  the  sensors.  A  highly  efficient  metamaterial 
structure  was  developed  to  provide  near  100%  absorption  at  3.8  THz  while 
simultaneously  serving  as  a  structural  layer  and  providing  access  for  external  optical 
readout.  The  fabricated  sensors  showed  responsivity  values  as  high  as  1.2  deg/pW  and 
time  constants  as  low  as  200  ms,  depending  on  the  configuration.  Minimum  detectable 

power  on  the  order  of  10  nW  was  observed,  demonstrating  that  the  sensors  can  operate 
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with  low-power  THz  sources.  Although  the  sensors  were  not  optimized  for  imaging,  the 
use  of  an  external  optical  readout  allowed  us  to  obtain  raw  images  of  the  QCL  beam 
indicating  the  potential  of  these  detectors  to  be  further  optimized  for  use  in  focal  plane 
arrays  for  real  time  THz  imaging.  The  successful  design,  fabrication,  and  characterization 
of  bi-material  MEMS  sensors,  using  metamaterial  absorbers  demonstrates  the  potential  of 
this  approach  for  THz  imaging. 

The  final  phase  of  the  research,  fabrication  of  THz-absorbing  Al/SiOx/Al 
metamaterial  layers  with  a  thin  Ti  conductive  layer  inserted  into  the  dielectric  spacer  to 
be  used  as  bolometric  layer  was  explored.  The  introduction  of  Ti  reduced  the  magnitude 
of  the  absorption  peak  when  placed  near  the  patterned  A1  layer  in  addition  to  a  small 
redshift.  The  measured  absorption  spectra  are  in  reasonably  good  agreement  with  those  of 
the  simulations.  The  results  indicate  that  the  bolometric  layer  should  be  placed  near  the 
ground  plane  for  achieving  maximum  absorption.  FE  simulations  of  the  effects  of  the 
bolometric  layer  on  absorption  suggest  weak  position  sensitivity  when  the  sheet 
resistance  of  the  bolometric  layer  exceeds  about  1  kQ/sq.  The  findings  indicate  that 
metamaterial  THz  absorbers  should  be  able  to  improve  the  perfonnance  of 
microbolometer  cameras  in  the  THz  range  without  compromising  their  thermal  time 
constants. 

A.  FUTURE  WORK 

The  knowledge  gained  in  this  work  can  prove  useful  for  the  fabrication  of 
microbolometer  pixels  with  integrated  THz  metamaterial  absorbers.  Therefore,  the  next 
step  in  this  research  should  be  the  fabrication  and  characterization  of  microbolometer 
pixels  using  the  optimized  metamaterials  with  embedded  bolometric  layers  for  THz 
imaging.  To  detect  a  change  in  resistance  in  a  microbolometer  pixel,  a  Wheatstone  bridge 
is  employed  in  readout  systems  [10]  as  illustrated  in  Figure  33.  To  ensure  the  load 
resistors  are  of  similar  resistance,  they  can  simply  be  other  microbolometer  pixels  that  are 
at  ambient  temperature.  Self-heating  is  an  inevitable  issue  with  a  microbolometer 
readout,  given  that  any  power  dissipated  while  probing  the  Wheatstone  bridge  will 
change  the  resistance  of  the  microbolometer  pixels.  Therefore,  a  possible  configuration 
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illustrated  in  Figure  33  has  two  thermally  shorted  pixels,  one  pixel  that  is  blind  to  THz  by 
joining  the  squares  together  (but  maintains  the  same  thermal  capacitance  as  the  sensitive 
pixel). 


Figure  33.  Wheatstone  bridge  configuration  for  THz  metamaterial  microbolometer 
characterization.  The  thermally  shorted  pixels  are  not  released  from  the 
substrate  and  the  blind  pixel  has  a  single  metal  square  instead  of  the  patterned 
layer  and  thus  will  not  absorb  THz  efficiently,  but  will  have  essentially  the 
same  thermal  properties  as  the  sensitive  pixel.  Several  Wheatstone  bridges 
are  fabricated  together,  their  common  ground  is  not  visible. 


This  configuration  will  automatically  cancel  much  of  the  self-heating,  allowing 
for  an  easier  readout.  Additionally,  a  separate  configuration  with  three  thermally  shorted 
pixels  and  one  sensitive  pixel  can  be  used  to  evaluate  the  self-heating  during  readout.  In 
commercial  systems,  an  array  can  be  read  out  pixel  by  pixel  with  the  same  Wheatstone 
bridge,  but  such  readout  systems  require  more  fabrication  steps  to  incorporate  the 
necessary  electronics. 
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